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Abstract. We consider a discrete polynuclear growth (PNG) process and 
prove a functional limit theorem for its convergence to the Airy process. This 
generalizes previous results by Prahofer and Spohn. The result enables us 
to express the Fi GOE Tracy- Widom distribution in terms of the Airy pro- 
cess. We also show some results, and give a conjecture, about the transversal 
fluctuations in a point to line last passage percolation problem. 



1. Introduction and results 

1.1. Discrete polynuclear growth. Recently there has been interesting devel- 
opments concerning certain special 1 + 1 dimensional local random growth models. 
This development has its starting point in the new results on the longest increasing 
subsequence in a random permutation, ||. We will not review all these devel- 
opments here. In this paper we consider a certain discrete growth model called 
the discrete polynuclear growth (PNG) model, p2| , a special version of which is 
closely related to the last-passage percolation problem studied in fll5| . It is a dis- 



crete version of the PNG model studied by Prahofer and Spohn, |29 , which can 



be obtained as a special limiting case. In the paper we will extend the results in 



29 1 to the present model and prove a stronger convergence result. We also obtain 
some preliminary results on the transversal fluctuations in the point to line version 
of the last-passage percolation problem, which should have many similarities with 
the corresponding problems for first-passage percolation and directed polmers. 

The discrete polynuclear growth (PNG) model is a local random growth model 
defined by 

(1.1) h(x, t + 1) = m&x(h(x - 1, t),h(x, t), h(x + 1, t)) + u(x, t + 1), 

x £ Z, t S N, h(x, 0) = 0, x £ Z. Here oj(x.t), (x.t) £ Z x N, are independent 
random variables, see ]22| . Typically they could be Bernoulli random variables. 
We should think of h(x,t) as the height above x at time t, so x — > h(x,t) gives an 
interface developing in time. We will treat a special case where lo(x, i) = if t — x 
is even or if |x| > t, and 

(1.2) w(i,j)=uj(i-j,i + j-l), 

£ Z?j_, are independent geometric random variables with parameter Oj&j, 

(1.3) = m] = (1 - a^)M,)"\ 

m > 0. We will mainly consider the case when Oj = b L = Jq, < q < 1, i > 1, and 
we we do this in the rest of this section. If we define 



(1.4) 



G(i,j) = h(i-j,i+j - 1), 
i 
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E it follows from ( [□]) that 

(1.5) G(t,i) = max(G(i - 1, j), G(i, j - 1)) + w(i,j), 

see proposition 3.10| . This leads immediately to a different formula for G(M,N), 

G(M,N) =max V* u;(i,j), 
(i,j)e7r 

where the maximum is taken over all up/right paths from (1,1) to (M,N). We 
can think of G(M, N) as a point to point last-passage time. It is also natural, from 
the point of view of directed polymers for example, to consider the point to line 
last-passage time, 

(1.6) G p i(N) = max G{N + K,N — K). 

\K\<N 

This makes it r easo nable to study the process K — > G(N + K,N — K), —N < K < 
N, which, by flO]), is the same as -fT — ► h(2K,2N — 1), i.e. the height curve at 
even sites at time 2N — 1. 

Let .Fi and F 2 denote the GOE respectively GUE Tracy- Widom largest eigen- 
value distributions, (3^]. It is known, JIB] , that there are constants a = 2y/q(l — 
y/q)' 1 and d given by (Q below, such that F[G(N,N) < aN + dN 1 / 3 ^ -> F 2 (0 
as N -> oo, and, fl, P[G pi (iV) < aiV + diV 1 / 3 ^] Fi(^) as JV ^ oo. Also, if 



the maximum in (1.6) is assumed at some point Kjy, which need not be unique, 
we expect Kn to be of order iV 2 / 3 , i.e. the transversal fluctuations are of order 
TV 2 / 3 . This can be seen heuristically, |22| , and there are some rigorous results for a 
related question, jl6) , j| , . These scales motivates the introduction of a rescaled 
process t — * Hiy(t), t 6 E, defined by 

(1.7) G(7V + «, iV — it) = 2v ^ iV + dN 1/3 H N ( l ~ ^ duN- 2/3 ), 

i - V? 1 + \/« 

and linear interpolation, |u| < iV, compare with [ p9[ . This is our rescaled discrete 
PNG process. The constant d is given by 

1 - 9 

In the limit when g is small and N is large, we can obtain the continuous PNG 
process studied by Prahofer and Spohn, j2{|. We want to extend their results to 
the present discrete setting and also prove a stronger form of convergence to the 
limiting process, a functional limit theorem. Before we can state the theorem we 
must define the limiting process which is the Airy process introduced by Prahofer 
and Spohn, @. 

We will approach by considering it as the top curve in a multilayer PNG 
process, compare and ^9|. This will lead to measures of the form introduced in 
sect. 1.2 and we will be able to use the formulas for the correlation functions derived 
there. The same methods can also be applied to Dyson's Brownian motion, compare 
with |l3| ], which can be obtained from TV non-intersecting Brownian motions. The 
appropriately rescaled limit as N — > oo of the top path in Dyson's Brownian motion 
converges to the Airy process, see below. This gives some intuition about what it 
looks like. Its precise definition is more technical. 



(1.9) A(t,£;t',0 
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The extended Airy kernel, p4| , [ pg| , is defined by 

f / °° e- A ^-^')Ai (f + A)Ai (£' + A)dA, if r > r' 
(too e~ A(T_T,) Ai (£ + A)Ai {g + X)dX, if r < r'. 

where Ai (•) is the Airy function. When r = r' the extended Airy kernel reduces 
to the ordinary Airy kernel, p3| . 

We dehne the Airy process t — > by giving its finite-dimensional distribu- 
tion functions. Given £1, . . . , £ m S K and ri < • ■ ■ < r m in R we define / on 
{n, . . . ,r m } x K by 

/(rj,a;) = -Xfe,oo)(a;)- 

Then, 

(1.10) P[A(n) < Ci, • ■ ., A(r m ) < £m] = det(/ + / 1/2 A/ 1 / 2 ) L2({rii ...^ }xR) , 

where we have counting measure on {ri,...,r m } and Lebesgue measure on K. 
The Frcdholm determinant can be defined via its Fredholm expansion, see sect. 
2.1 below. We will prove in section 2.2 that f 1 / 2 Af 1 / 2 is a trace class operator 
on L 2 ({ti, . . . ,r m } x R), so this is also a Fredholm determinant in the sense of 
determinants for trace class operators. Note that in particular 

(1.11) P[A(r)<e]=i f 2(0- 

This defines the Airy process. It is proved in |^9f that it has a version with contin- 
uous paths, which also follows from the results below. As mentioned above, another 
way of understanding the Airy process is as follows. Let X(t) = (Ai(t), . . . , Ajv(i)) 
with Xi(t) < ■ ■ ■ < Ajy(i), be the eigenvalues in Dyson's Brownian motion model, 
||, for GUE with stationary distribution Z^A^A) 2 ]jj =1 exp(-A 2 ). Then, 

lim V2N 1/e (X N (N~ 1/3 t) - V2N) = A(t) 

JV— ►co 

say in the sense of convergense of finite-dimensional distributions. This can be 
proved using the methods of the present paper, and using techniques from | jr9[ , it 
is possible to get an integral formula for the (extended) correlation kernel. The 
details will not be given here. This scaling limit has been studied before, see [ fl3| 
and references therein. 

In analogy with the results of |2^] , we can show that the rescaled height process 
Hn converges in finite dimensional distributions to the Airy process. 



Theorem 1.1. LetH^ be the process defined by (1.7). Then for any fixed ti, ... ,t m 
and £i, . . . ,£ m , 

(1.12) lim F[H N (h) < £i, ■ ■ -,H N (t m ) < U] 

P[A(ti) < £i + A, ■ ■ ■ , A{t m ) <U + t 2 J, 
where A is the Airy process. 

This result can be sharpened to a functional limit theorem. 
Theorem 1.2. Let A(t) be the Airy process def ined by its finite- dimensional dis- 



tributions, (L1C). Also, let -ffjv(^) be defined by (LV) and linear interpolation. Fix 
T > arbitrary. There is a continuous version of A(t) and 

(1.13) H N (t)^A(t)-t 2 , 

as N — > oo in the weak* -topology of probability measures on C(— T, T). 
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The theorem will be proved in section 5.2. As a corollary to this theorem and 
the results of Baik and Rains, [Q, we obtain the following result which expresses 
the GOE largest eigenvalue distribution Fi in terms of the Airy process. 

Corollary 1.3. For all £ e R. 

(1.14) F 1 (0 = P[sMA(t)-t 2 ) <e 



The proof of (1.14) is very indirect. It would be interesting to see a more straight- 
forward approach. 

As discussed above we are also interested in the transversal fluctuations of the 
endpoint of a maximal path in the point to line case. In our discrete model this 
is not well-defined, there could be several maximal paths. Consider the random 
variable 

(1.15) Kn = inf{u ; supi?Ar(£) = sup -£/at(£)}, 

t<u tea 

the first point that gives the maximum. The corresponding quantity for the limiting 
process H(t) = A(t) — t 2 is 

(1.16) K = mi{u; supH{t) = sup H{t)}. 

t<u teif 

We would like to show that Kn converges to K so that we could call the law of 
K the asymptotic law of transversal fluctuations. Unfortunately we can only prove 
this under a very plausible assumption on the Airy process. We can show, 

Proposition 1.4. The sequence of random variables {Kn}n>i is tight, i.e. given 
e > there is a T > and an Nq such that 

F[\K N \ > T] < e 

for all N > N . 

The assumption we need to make on the Airy process can be formulated as 
follows. 

Conjecture 1.5. Let H(t) = A(t) — t 2 . Then, for each T > 0, H(t) has a unique 
point of maximum in [-T, T] almost surely. 

If we accept this we can prove 

Theorem 1.6. Assume that conjecture \l.5\ is true. Then Kn —> K in distribution 
as N — > oo 

The law of K is thus a natural candidate for the law of the transversal fluctua- 
tions. It would be interesting to find a different, more explicit, formula for this law. 
Assuming the truth of the same conjecture it may also be possible to prove that 
the endpoints of all maximal paths, or asymptotically maximal paths, converge to 
the same limit K . 

By using the limit results of Q, proposition 3.12 and theorem [3.14 we can ob- 



tain the correlation functions of the eigenvalues of the succesive minors 
{hij)i<i.j<ki 1 < k < N, of an N x N GUE matrix H = Qi%j)i<i^<.N • ln this wa Y 
it is possible to get the Airy process as an appropriate limit of the succesive largest 
eigenvalues of . More details will be given in future work. 

We could also get the Airy process by looking at the largest eigenvalues of coupled 
GUE-matrices, which is similar to looking at Dyson's Brownian motion model for 
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GUE. In [£9| Prahofer and Sp ohn r aised the problem of finding differential equ ation s 
for probabilities of the form ( 1.1C ) generalizing the Painleve II formulas for ( 1.11 ). 
In jjj the spectrum of coupled random matrices is studied and it would be interesting 
to see if the results of this paper shed some light on this problem. 



1.2. Measures denned by products of determinants. Probability measures 
given by products of determinants has been studied in several papers, e.g. by 
Eynard and Mehta, j[o|, in connection with eigenvalue correlations in chains of 
matrices, by Forrester, Nagao and Honner, fll3| , in connection with Dyson's Brow- 
nian motion model and by Okounkov and Reshetikhin, |2^| , when introducing the 
so called Schur process. The problem is to compute the correlation functions and to 
show that these are given by determinants so that we obtain a determinantal point 
process, The same type of correlation functions are also obtained by Prahofer 
and Spohn, [^9), in a cascade of continuous polynuclear growth (PNG) models. We 
will study a class of measures which include all the above as special cases and show 
that we obtain determinantal correlation functions. As an example of the result 
we will in sect. 2.3 investigate random walks on the discrete circle using the same 
strategy. This will lead to an extended discrete sine kernel, compare with [^8). We 
will see in sect. 3 that our main topic the discrete PNG problem fits nicely into 
this framework. This particular application is very close to the Schur process in 
and their results could also have been used. In fact, we rcdcrivc their main 



formulas. 

For r G Z let x r = {x\, ...,<) G R" and x = ( X - M +\ x M ^), M > 1. We 
think of a; as a point configuration in {— M + l,...,M-l}xK, and we also specify 
fixed initial x~ M and final x M positions. Let 4> r .r+i : I 2 ^ C, r 6 Z, be given 
transition weights. The weight of the configuration x is then 

M-l 

(1.17) w nM {x)= [] det(<^, r+1 «,4 +1 ))^ =1 . 

r=-M 

Let dfi be a given reference measure on R, typically Lebesgue measure or counting 
measure. We assume that \^ r +i(x,y)\ < c r (x)d r (y), where c r G L 1 (M, /i) and 
d r G L°°(R, /i), —M < r < M. This assumption is not necessary but is convenient 
and suffices for the convergence of all the objects we will encounter. The partition 
function is 

(1-18) Z nt M = , ,J M ; / w nM (x)dn(x), 

where d/i(x) = Ylr=—M+i TYj=i ^(^P- ^ e w ^ n assume that Z n _M ^ so that we 
can define the normalized weight 

(1-19) Pu,m(x) = . ,, 2 m_ 17 w n>M (x). 

(Til) ^n,M 

If w Ut M(x) > 0, this is a probability density on (R n ) 2M_1 with respect to the 
reference measure d/J,(x). The (k-M+i, ■ ■ ■ , fc/y-i)- correlation function can now be 
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defined in a standard way by 

( i 201 Ri u (t~ m+1 t - m +! r M -i T M -i ) 

yi.tu) Jik_ M+1 ,...,k M - 1 (^i , ■ ■ ■ i^k-M+V ' ' " ' 1 k M -i> 



-M+l ' 

where fc = fc-ju+i + • • • + fcM-l> < fc., ■ < n. 

Given two transition functions we define their convolution by 

4>*ijj(x,y) = I (j>{x,z)il>{z,y)dn{z). 
Jm. 

Set 

(f) r<s (x,y) = Oryr+l * ■ • • * (f) s -l iS ){x,y) 

if r < s and r)S = if r > s. Let ^4 = (Ay) be the n x n matrix with elements 
Aij = 4>-m m{xJ m , xf 1 ), 1 < i,j < n- By repeated use of the Heine identity: 
(1-21) 

^ f detOifo-))"^! det^iixj^^d^x) = det ( f &(f)^-(t)d/i(t) ) 

we see that Z„_m = det A. Hence det A ^ by our assumption. Define a kernel 
K n,M . (|_ M + i i ... ] M-l})xl) 2 ^Cby 

(1.22) K n > M (r, x- s, y) = K n ' M (r, x; s, y) - <j>r iS (x, y), 
where 

n 

(1.23) K n ' M (r,x;s,y) = £ r , M (z, ){A- 1 ) ij tl > - Ml .(xJ M , y). 

In the case M = 1 the kernel X has appeared before, see @, and also @. 
Theorem 1.7. TTie correlation functions defined by (1.200 are 9^ ven by 

^1. nk_ u+u ...,k M -i\- L l ■• ■ ■ ■ i x fc_ M + i i ■ • • i x l J ■ ■ • > x k M -i' 

= det(K n - M (r, «L; s » a! i,))-M<r,s<Af,o<tr<kr,o<i«<k a 



The determinant in the right hand side of (1.24) has a block structure with the 
blocks given by r, s and having size k r x k s . The theorem will be proved in section 
2.1. 

A case of particular interest is when the transition weights are given by Fourier 
coefficients. We are then in a situation similar to that in Let f r {e 10 ) be a 

function in L : (T) with Fourier coefficients f r . Assume that the transition weights 
are given by 

(1-25) <f> r>r+ i(x,y) = f r (y- x), 

—M < r < M, x, y € Z and that the initial and final configurations are given by 
Xj M = xj 1 = 1 — j, j = 1, . . . , ri. If we set 

s-1 

fr,s( Z ) = II -feW' 

z = e , then 

(1-26) (j> r , s (x,y) =f r>s (y-x) 
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for r < s. The matrix A defined above is then a Toeplitz matrix with symbol 

M-l 

(1-27) a(z) = f- MtM {z) = [] 



--M 



Define 



(1.28) K%?(z,w)= £ K n > M (r,x;s,y)z*w-y, 



where K n ' M is given by (1.23). When the transition functions and the initial and 
final configurations are given in this way we are able to give a formula for the limit 
of this generating function as n — > oo. 

Proposition 1.8. Assume that f r (z) has winding number zero, a Wiener-Hopf 
factorization f r (z) = /+(z)/~(z) and is analytic in 1 — e r < \z\ < 1 + e r for some 
e r > 0. Furthermore, suppose that 

^2 M a |«nl < °°i 

nGZ 

for some a > 0, where a n are the Fourier coefficients of the symbol a{z) given by 
(1.27). Set e = mine r and 

(1.29) iC™(z,w) = ^—G(z,w), 



z — w 



where 

(1.30) G(z,w) = 



T\t=r ft C z ) Tlt=-M ft~(w) 

rit=-M rit=s ft (~t,) 



Then, for 1 - e < \w\ < 1 < \z\ < 1 + e, 
(1.31) 



Furthermore, 



< \fr,M(j)\\f-M,s(^)\ (_L + \ w \n/2 



(|«|-1)(1-H) \n a 1 1 \z\ n / 2 



(1.32) ^(x, V) = ^f e^-^ B G{e ie ,e ie )d6, 

for r < s. 

The same type of formula for the limiting kernel was obtained in f28|| . The 
formula will be proved in section 2.1. This proposition makes it possible to compute 
the asymptotics of the kernel given by ( 1.22f ) in certain cases, since it gives an 
integral formula for the n — > oo limit of K n ' M . 

The outline of the paper is as follows. In sect. 2 we will give some general results 
for measures of the form (1.17) and then as an example discuss nonintersecting 
random walks on the discrete circle. The next section applies the general theory to 
the discrete PNG model and gives more explicit formulas. In sext. 4 asymptotic 
results for the correlation kernel appearing in the PNG model are stated and proved, 
and these are then applied in sect. 5 to prove the necessary estimates needed for 
the functional limit theorem and the transversal fluctuations. 
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2. Determinantal MEASURES 

2.1. General theory. In this section we will prove the results of section 1.2. We 
will prove theorem [O] using a generalization of the method of |34j , for (3 = 2 
random matrix ensembles, see also ]lq] . It is also possible to generalize the approach 
of |^0|, which is closer to the original Dyson approach. 

Let Am = {—M + 1, . . . , M — 1} x R, A be the counting measure on {— M + 
1, . . . , M — 1} and v = X ® ^i. Furthermore, we let g : Am — * C be a bounded 
function and define 

If n 
Zn,Af\g] = , n2 M-i / IT Y[0- + 9{r,x r j ))w niM (x)dfj,(x). 



We want to compute Z n ^M\g\l Z n ,M$\- Using the Heine identity (1.21) repeatedly 
we see that 



1 



,m[q] = i .S2M-1 / det(^_ M ,-M+i(a; i ' ,x ' ))i< lJ <„ 

M-l 

J| det((l - g(r,x[))0,. !r+ i(a;[,a;^ +1 ))i<i J < n d^(5) 



M-l 

X 

r=— M+l 



det 



/ t-M.-M+lfa M ,t-M+l) TT ( 1 +9(r,t r )) 

y^'- 1 k | <M 



Af-2 



X 



Y[ </>r,r+l(U,U+l) <j>M-l,M < \ t M-\,xf)d lM 1 jJ,{t) 



l<ij'<n 



Now, 



2 M-l 

J| (1 + 9{r, U)) = 1 + Y 9{ri,t ri ) ■ ■ -g(re,t re ), 

\r\<M (=1 -M<ri<—<rt<M 

and hence 

2M-1 

-M 



/ 2M-1 

(2.1) Z„ >M [ 9 ]=det X] Yl / 

\ £=1 -Af<n<— <r7<Af' H< 



-M,rx{%i ,tl) 



Y\_9i r s,ts)4 l r a .r s+1 ( t s,ts+l) 9(n,h)(f) re ,M{tl, xf)d l fJ,(t) , 
\s=l / I \<i,j<n 

where we have used the notation of sect. 1.2. If we set g = we obtain Z Mi m[0] = 
Z n ,M — det A as before. B y de finition <j) r ^ s = if r > s, and hence we can remove 



the ordering of the r^'s in (2.1). We find, 



2M-1 



(2.2) 



det Uj +Y,(A- 1 ) ik ]T ]T I <t>-M, n {xl M ,t x ) 

\ i, — i 6 — i nrf"^-*. ^ns"^ 



z n ,M[g] 

Zn ' M ^ \ fc=l £=1 — M<r m <M * 

|nj( r " t «)^.,fa +1 ( t »i W)^ 9(n,k)4>r e ,M( t i,Xj I )d i fj,(t) 



l<i,j <n 
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Write ijj(u, t;v,s) = g(u, t)4> u ,v(t, s), and define ip*°{u, t;v,s) = 5 uv 5{t— s), ip* 1 = ip 
and 

?A* (r+1) (u, t; v, s) = / ip(u, t; mi, ?i)^(mi, £ x ; m 2 , &) • ■ ■ <K"V, £r! u, s)d r v(m, f) 

for r > 1. Note that, since <£ ns = if r > s, we have V* € = if t > 2M - 1. This 
follows immediately from the definition. The formula (2.2) can now be written 

(2.3) 

Zn,M[g] 



= det I S i:j +y2(A- 1 ) ik [ dv{uX) I dv{v,ri) 

V fe= l ^ ^A m 

/2M-1 \ 



If K(x.y) is an integral kernel on L 2 (f2,^) we define the determinant det(I + 
K)l 2 (q.,il) via a Fredholm expansion, 

°° 1 f 

(2.4) det(J + lO L 2 (niM) = XI — / det(if(a; l ,2: J )) 4J=1 ,... !m d m M(x). 

We assume that if is such that all the integrals are well-defined and the se- 
ries converges. For example, by Hadamard's inequality, it is sufficient to require 
that |-K"(a;,y)| < a(x)b(y), where a e ^(fy/x), b G L°°(0,/i). Note that if 
f2 = {1, . . . , n} and fj, is counting measure this is the ordinary determinant det(5ij + 
K(h j))i,j=i, Let K(x,y) be an integral kernel from L 2 (f2i,/ii) to L 2 (f22,A*2) 
and L(x,y) an integral kernel from L 2 (H2 , ^2) to L 2 (i!i,/ii). Then 

L * K(x, y) = I L(x,z)K(z,y)dn 2 {.z) 
Jn 2 

is an integral kernel on L 2 (Oi, /hi). Furthermore, 

(2.5) det(J + L * AV(n llMl ) = det(/ + if * L) L 2 (naiMa ). 



This is easy to see using the Heine identity in the definition (2.4). 
Set 



fc=i 

- 2M-1 



so that, by (|2j) and (g. 



det(<5y + (6*c)(i, j))i<i,j< n 



= det(I + c*b) L 2 {AM v) . 
Now, a computation shows that 

2M-1 
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where K is defined by ( 1.23 ). Thus, 



(2.6) 



,M 



2M-1 



[0] 



det(7+( J2 r ie - iy )*(9K)) L : 



(Am.i')- 



1=1 



The kernel in ( |2.6| ) has finite-rank so the sum (2.4) in the definition of the deter- 
minant actually has finitely many terms. We now claim that the right hand side of 

(2.6) equals 

det(J- tl> + gK)L 2 (A M ,u), 
which is what we want. Formally the computation goes as follows. The expression 
in (2.6) is det(J — (I — ip)~ 1 gK) and we multiply this by det(J — tp) — I. Since we 
are only working with determinants defined by a Fredholm expansion the product 
rule is not obvious, so we will give a proof in this special case. 
Write a = gK. We will prove that for any z, w 6 C, 
m 

(2.7) det(I + w^2z : >ip* i3 - 1] *a) L 2 (AjVf) = det(/ - zip + zwo) l ^ m) , 

j=i 

where m = 2M — 1. The left hand side is a polynomial in z,w so it suffices 
to prove (2/7) for sufficiently small. In that case, under our assumption 

on the (f> rtr +i, all the expressions beloware well-defined and convergent. Write 
b = J2T=i * a. Then, see e.g. ||, 



det(7 -I- wb) L 2 



(A,, 



ex p(E 



(-1)^+1^ 



fc=l 



b* k {t,t)dv{t)) 



Am 



and 



dct(7 + z(-ip + w))f,2( Ajf) = exp(J^ 



(-l) k+1 z k 



k=l 



i) + wo)* k (t,t)dv{t)). 



Set c = za, d = zip. It suffices to show that 



E 



k=l 



(2.8) 

The equality 



E 

k=l 



k 

(_!)fc+i 



(t,t)dv(t) 



{~d + wc)* k {t,t)dv(t). 



Am 



holds for w — since 

{-dy k {t,t)dv{t) = {- z ) k f ip* k {t,t)dv{t) = o 

if k > 1. This follows from r?s = for r > s. Hence it is enough to show that the 
derivatives of the two sides of (2.8) coincide, 

m— 1 



oo „ I m— 1 

E 

fe=0 JAm \j=0 



*(fe+l) 



(i, *)di/(t) 



OO „ 

53 (-1)" / ((-d + u>c)*"*c)(t, *)<&/(*)■ 

n=0 ^ Aj w 
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To prove this last equality is a straightforward but somewhat tedious computation, 
which is based on expanding both sides and showing that the coefficient of w k is 
the same on both sides. We omit the details. 
We have proved 

Proposition 2.1. Let p n ,M(x) be defined by 3 j and assume that (j> rtT j r \(x,y) 
satisfies \<j> r>T+ x(x,y)\ < c(x)d(y) with c e L 1 (R, d e L°°(R,[i), -M < r < M. 
Furthermore, let Am = {—M + 1, . . . , M — 1} X R, v — A <X> \i, where A is counting 
measure on {— M + 1, . . . , M — 1} and let g : Am — > C be a bounded function. Then 



(2.9) 



l} 2 



J] Y[(l+g(x»))p n , M (x)dn{x) = det(I + gK) L , 



(A M , v) j 



where K is given by (1.2i), and the determinant is defined by using the Fredholm 
expansion (2.4)- 



Theorem |l.7|is a direct consequence of (J2.9D , compare with the discussion in J34{ . 

If AT m = {1, . . . , to} and A is counting measure on X m , then L 2 (X m , A) = R m and 
we have a chain of isomorphisms L 2 (X m x f2, A<8>//) = L 2 (X m , A)(g>L 2 (r2, /i) = R 5 * (g> 
L 2 (fl, jj) = L 2 (fl, //) © • • • ©L 2 (il, fi), where we have to terms in the last direct sum. 
We can think of an element in R > © L 2 (il, /i) as a column vector (fi(x) . . . f m (x)Y, 
where fi e i 2 (fi, /x), 1 < i < to. Hence, an operator on L 2 (AT m xfi,A®/i) defined 
by an intgral kernel K(r,£;r' , £') can be thought of as a block operator on these 
column vectors with block kernel {K(r, £; r', £'))i<r,r'<m.' 

We also want to prove Proposition [T^. Let us write T n (a) for the n x n Toeplitz 
matrix with symbol a and T(a) for the one-sided infinite Toeplitz matrix with 
symbol a. Consider the function K n - M (z,w) defined by (1.28) and let the symbol 
o be given by fll.27|). Then, 



(2.10) 



1 1 " 

-frM-)f-M, s (-) E ^Pr^v 



i>J=l 



To proceed we need a formula for the inverse of a Toeplitz matrix. We will use the 
following result which follows from theorem 1.15 and theorem 2.15, together with 
its proof, in ||. 

Proposition 2.2. Assume that a(z) = a+(z)a,-(z), z <E T, where 

oo oo 

(2.11) a+(z) a-(z) = E a =^~"' 



71 = 



71 = 
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S^Lo(l a n I + l a -rJ) < °°; an d that (a(z) has winding number zero. Furthermore, 
suppose that 



(2.12) 



H a \a n \ < oo 



for some a > 0, where a n is the Fourier coefficient of a{z). Using ( 2.1 f[ ) we can 
extend a + (z) to \z\ < 1 and a_(z) to {\z\ > 1} U {oo} and we assume that they 
have no zeros in these regions. Then, T n (a) is invertible for n sufficiently large and 
there is a constant C ( which depends on a) such that 



(2.13) \[T-\a)] jk - [T(a; 1 )T(aI 1 )] jfc | < Cmin( 



for 1 < j, k < n. 



1 



1 



(n + 1 - k) a ' (n+l-j) 



-) 



We can now give the proof of proposition 1.8. 



Proof, (of Proposition 1.8). The function a defined by (1.27) has a Wiener-Hopf 
factorization a — a + a_ where 



M-l 



t=-M 



and all the assumptions of the previous theorem are satisfied. By (2.13) 

n n 

z-^T-^a)]^ ^[Tia+MaZ 1 )}^ 

n 

< C £ M~>K min( — — — ) 

n + 1 — i) a (n + 1 — i) a 



< 



C 



— (— + \ w \ n / 2 + 1 ) 

(H-iKi-Mr^ m™ /2 



Also, 



^ z _, [T(ffl^ 1 )r(a: 1 )] 1J w J 



i > n or j > n 



< c 



\w\ n + \l/z\ n 

(N-i)(i-H)" 



Set 6± = l/a± and note that (6+)fc = if k < and (&-)& = if k > 0. We can 
now compute 



£ z- 4 [T(a; 1 )T(a: 1 )] 4 



oo 



fe=l \i6Z / 

1 w/z 
a + (l / z)a-(l / w) 1 — iw/z 
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It follows that 

7 1 1 - 

-frM-)f-M, s {-) E ^m^maz 1 )] 

in -7. in — ' 



ijW 3 

z to' ' — ' 



rit=r ft~d)ft (z) il=-m ft~(w)ft (w) 



rAf-l f+ ,l, ... ! ■ = «0, 

U=-Af/t Ij, 

and the proposition is proved. □ 
We have 

1 f dz f z 

(2.14) K(r,x;s,y) = K(r,x;s,y) = — — - / — - / w y - l dw G{z,w). 

(2my J lr2 z x+L J Jri z-w 

if r > s, where 1 — e < r\ < r2 < 1 + e. Using the residue theorem it follows that 
for r < s, 

(2.15) K(r, x; s,y) = ^ ^d W -^-G { z, W ). 
1 — e < n < r2 < 1 + e, compare p8[. 



2.2. The extended Airy kernel. The extended Airy kernel is defined by (1.9). 
We can also define a modification by 

/>oo 

(2.16) A(t^-t\0= e~ A(r ~ r ' ) Ai(C + A)Ai(e' + A)dA, 



which is well-defined both for r > r' and for t < t' by the following standard 
estimate for the Airy function, 

(2.17) |Ai(£)|<C M e- 2| « |3/2/3 

for £ > — M. Both A and A have a useful double integral formula. 



Proposition 2.3. The extended Airy kernel (l.i) is given by 



i r r e i£z+i£'w+i(z 3 +w 3 )/3 

(2.18) A {T^T>,i') = -— 2 dz dw T ,_ T + i{z + w) , 

where rj, rj' > and rj + fj 1 + t — t 1 < in case r' > r. ^47so ; t/ie modified kernel A, 
(2. It ), is given by the same formula but where we now require that rj+r/' +t—t' > 0. 

Proof. This is straightforward using the identities 

e -A(r-T -iz-iw) dx = 

r — t + i[z + to) 
if r — r' + n + rf > and 

1 



e 



t' — t + i(z + w) 

if t - t' + rj + -q' < 0. □ 



If we move the contour of integration between the two cases in proposition 2.3 
we pick up a contribution from the singularity and we obtain 

(2.19) A(t, £; t', O = A(r, £ r', £') - <£ T , T ,(£, £'), 
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where d> TT i = if r > t' and 



(2.20) 4> T y{^) = — c 
^47t(t' - r) 

if r < t'. Combining (p|), ( |2.16D and Q2.19| ) we see that 



-(?-e') 2 /4(r'-r)-(r'-r)(C+5')/2+(r'-r) 3 /12 



(2.21) 



4>r,r> & £') = / e- A ^-^') Ai (£ + A)Ai (£' + A)dA 



if r < t'. We would also like to show that the operator in ( 1.1C ) is actually a trace 
class operator. 

Proposition 2.4. Let f(r,x) bew a non-negative function in L°°(R) for each t G 
{ti, . . . , r m }, where tj. < ■ • • < T m . Assume also that /(r^, x) = if x < Mj~ for 
some number Mk, k = 1, . . . , m. Then, the kernel 

f^x^A^x-r^x^fir^x') 1 ' 2 

defines a trace class operator on L 2 ({t\, . . . , r m } x R). where we have counting 
measure A on {t\, . . . ,t,„} and Lebesgue measure fi on K. 

Proof. We will prove the result by factoring into two Hilbert-Schmidt operators. 
Let H (t) = 1 if t < and H(t) = if t > 0. Set 

/OO 
e -V(r-r') Ai ^ + y ) Ai + 
-OO 

For i < j wc define 

Bij(r,x;T' ,x') = B(t,x;t' \x')5 TfTi 5 T ^ Tj 



so that 



S(t, x;t',x') = Bij(T,x;T',x'). 

l<_i<j<m 



Since, by ( |2.19| ) and ( |2.2lj ), A = A - B, it suffices to show that / 1 / 2 i/ 1 / 2 and 
f^l^Bijf 1 ! 2 , 1 < i < j < m, are trace class operators. 
Set 

a(r, z; a, y) = -L/(r, x) 1/2 Ai (z + y) e ^ (r ^ ) X[o,oo)(y) 



6(<7, y; r', z') = -^X[o,oc) (l/)Ai (a:' + yK^^iV, x'f' 2 



Then a and 6 are Hilbert-Schmidt kernels on L 2 (A m , A(g>/i), A m = {ti, . . . ,r m } : 



We have 



A m J A 
1 



< 



a(r, x; cr, y)| 2 d(A (g) /l«)(t, a;)d(A (8 ^)(cr, y) 

f(r, x)Ai (x + y) 2 x % oo){y)e- 2y{T - a) dxdyd\{T)d\{a) 

V / / dyAi(x + y) 2 e 2, ^ Tm - T1 \ 



A m JA 
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where M — min(Mi, . . . , M m ). Using ( 2.17 ) we see that the integral in the last 
expression is < oo. The proof that b is a Hilbert-Schmidt kernel is analogous. Now, 

o(r, x; a, y)b(a, y; t' ,x')d{\ ® (i)(t, y) 

- V f^xY^fiT^x') 1 ' 2 e-y^')Ai(x + y)M(x'+y)dy 
m * — ' ; n 

CT e{ri,...,-r m } J ° 

= f(T,x) 1 ^A(T,x;r\x')f(T',x / ) 1/2 - 

Hence, the operator f 1 / 2 Af 1 / 2 is trace class. 
Next, set 

c 13 (t,x;t',x') = -Lf(r,x) 1 l 2 ki{x + y)e-y^-^l 2 5 T .^ 



(it is independent of a) and 

di^x-y ,x') = ^/(r',x') 1/2 Ai(x' + y)e-^- T ')/ 2 ^, T3 . 



Then, 



c tj (T, x; t', x')dij(T, x;t' ,x')d{\ ® fi)(a, y) 

•I poo 

= - V f(r, xf^fir', x') 1 ' 2 / e-y<r-r >Ai (x + y)Ai (x' + y)dyS T , 1 ,)..- 
= Bij(T,x;T',x'). 

It remains to prove that cy and dij are Hilbert-Schmidt kernels. Consider c^; the 
proof for dij is similar. We get 



1 £ 

m z — ' 



o",r£{ri,...,T m } 

f OO /"DO 



oo /»oo 



/(r, x)Ai (i + y) 2 e- v(T - T ^dxdyS TtTi 



OO «/ — OO 



— oo «/ — OO 



< 



< 



/(Ti,aj)Ai (a: + y) 2 e- y(Ti - T ^dxdy 
dx I dyAi(x + y) 2 e y{T J- n) 



Mi J -oo 
poo 

dx / dyAi(x + y) 2 e y ^- n) 

M Jo 

x , da; / dj/Aifc + j/jV^-^. 



The first integral in the last expresion is < oo by ( 2.17 ). Now, by ( 2.17 ) 



Ai(x + y) 2 dx = Ai (xfdx < (7(1 + \y\) 

Mi J Mi+y 
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since the Airy function is bounded. Hence, 



r° 

dx / dyAi(x + yfe y( - T i- Ti) 

Mi J-oo 

<C||/|U / (l + |y|)e^-^dt/<oo ! 

since Tj — n > 0. This completes the proof. □ 

2.3. An example: random walks on the discrete circle. We will consider non- 
intersecting walks on the set Zn of integers modulo N, the discrete circle. This 
type of model has been analyzed in |l2] ] and we will show how it fits into the present 
formalism. We have 2M — 1 copies of Zjv, where the first and the last are identified 
so that we have periodic boundary conditions in the time direction. We will have 
are non-=intersecting paths on the discrete torus. Let x r S Z^ be the particle 
configuration (n particles) on the r:th discrete circle, \r\ < M,x~ M+1 = x M ~ x . 
Assume that n is odd, n = 2z/ + 1 and that the transition probabilities for the walks 
are given by 

!P, i£y-x = l 
q, i£y-x = 0. 
0. otherwise 

for x, y € Zjv, where p, q > and p + q = 1. The transition probability for non- 
intersecting paths from a configuration x r to a configuration x r+1 is 

det(^ rjr+ i(a:[,^ +1 )i< ij < n ). 

Write x = (x~ M+1 , x M ~ r ) G (Zjv) 2 *^ 1 for the total configuration. The prob- 
ability of x is 

M-2 

(2.22) q n ,N,M{%) — II det(</) rjr+ i(a;[,a;J +1 )i<j J < n . 

r=-M+l 

We will use discrete Fourier series on Zjy, 



where z — e 27ri f N . Also, we can represent Kronecker's delta on Zjv as 



(2-23) 5 xy = 1 



z 



k(x-y) 



Let Zjy- = {x e Zjy- ; < x\ < ■ ■ ■ < x n < N} be all ordered configurations of n 
particles on Z N . If ar -M+1 ,a: M ~ 1 G Z™ , then 

det^-M+l x M_1 )l<2j'<n = ^- M + 1 ,l M " 1 
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by the definition of the determinant. This determinant can be rewritten using 
(2.25) and Heine's identity, 



(2.24) 



+i jB «-i = det(- }2 z ' 3 } )i 



— ^ det(z fc ^ )i<i,j< n det( 
fci,...,fc„ez„ 

if x _M+1 ,:r M_1 G Zjy. This leads us to the measure 

1 



<2 j'<n 



Z ' 3 )l<i,i<r 



J N ■ 

Pn,N,M(%) 



-Qn,N,M(x)5 x -M + l x M-l 



(n!) 2M - 1 Z„, JV , M - 
where ^n,w,M is the normalization constant. 

Let g(r,x), \r\ < M , x € Zjv, be given functions and set 

n 

g(x)= n n( i+ ^^))- 

|r|<Mj=l 

We want to compute the expectation 



(2.25) 



G(x)p n>N> M(x) 

ee(z™)2M-i 



N n Z, 



n,N,M 



^ G(x)w n ,N,M{k;x), 



where 
(2.26) 



w n ,N,M{k;x) = aet{z 1 3 )q n ,N,M\ x ) det ( 2 * ' ) 

M-l 

= J] det(^, r+1 «,^ +1 )). 

r=-M 

Here we have set (j>-M,-M+i(h, X J M+1 ) — zh ' Xj i ^M-i,M(if _1 , %) = 2" 
and £~ M = xf^ = foj € Zat. We have a measure of the form ( p.. 17 ). Set 

(2.27) Zn t N !M (k) = ^ U>n,N,M(k;x) 



and note that G = 1 in (2.25) gives 



(2.28) 

Let us also write 



TL\ x — > 

Z n ,N,M = 2^ Z n ^ N ^ M (k). 



Pn,N,M(k,x) 



(nlf M ^Z, hNM (k) 
The expectation fl2.25| ) can then be written, using ( p. 28 ) 

Z n ,N,Ai(k) 



w n ,N,M{k\x). 



(2.29) 



a- e 



EfeGZ" Z n,N,Al(k) 



E ny N,M(k': G) 



where E n ^ t M(k; G) is the "expectation" of G with repect to the measure p n ,N,M{k, x) 
This "expectation" can be computed using the standard framework. Let f(n) be 
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equal to q if n = 0, p if n = 1 and otherwise, n G %n, so that (j) r . r +i(x,y) — 
f(y - x), for —M < r < M — 1. Then f(£) = q + pz l , I G Zjy, and by standard 
properties of convolution 

(2.30) <p r , s (x,y) = f s - r (y-x), 
—M < r < s < M — 1. From this we see that 

A l3 = <P- M<M (ki, kj) = ]T z kiX f 2M - 2 (y - x)z- k >y 

= Nf(N - h) 2M - 2 S k ^ k] = Nf(N - fc,) 2M ~%- 

if keZ^. Thus, 

(2.31) A = (AT/(]V - k t ) 2M ~ 2 6 l3 ) lJ=h ..., n . 
Now, 

(2.32) Z ntNtM (k) = detA = N n {](<? + p e - 2 ^' N ) 2M ~ 2 . 

i=i 

This is always non-zero if p ^ q. If p = q = 1/2, then we assume that N is odd, 
which also implies that det A ^ 0. We obtain 



(2.33) K n (k;r,x;s,y) = £ ( £ " ) ^ W - ki) 2 -™ 6^ 

i,j=l \£6Zjv / 

VmSZjv / 
1 ™ 



N 

i=l 



where we have indicated the dependence of the kernel on k. Note that this kernel 
is independent of M. We have 

1 " 

(2.34) K n (k; r, x;s,y)^-^2 f( N ~ h) s - r z kAy ~ x) ~ <Mx, y), 



N 

i=l 



where 



(2.35) ^(^)4E/r^ (rl 



TV 

fez 



if s > r an d <fi r , s (x, y) — if s < r. Computations similar to those leading up 
to formula ( 3.27 ) below show that if we assume that g(r, ■) = if \r\ > M q, then 
E n ,N.Ai(k,G) = E n ^ NiMo (k,G) for M > M . Hence, the expectation ( [2.29 ) can be 
written 



Elim 
\ M->oc 



r2.:VU > ^ ( lim -WW \ E {k>G y 

+0 ° LfceZ™ Z n.N.M(k) I 



Lemma 2.5. Let on — i — 1, i = 1, . . . , v + 1, a2i/+2-i = N — i, i = 1, . . . ,v, 
n = 2v + l. Ifke Z%, then 

to Q'7^ r Z n N M {k) 

(2.37) lim == 7TT=6a,fc. 
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Proof. We use the explicit formula fl2.32Q for Z ni N,M(k), 

n 

Z n , NM (k) = N n Y[(q + pe-™**/*)™-*. 
3=1 

Now, 

\q+pe- 2 ^/ N \ 2 = p 2 +q 2 + 2pqcos^^. 
This is maximal (= 1) if ki = 0(= N) and it is easy to see that 

(p 2 + o 2 + 2pgcos -^) < } [{p 2 + q 2 + 2pqcos —^-) 

3=1 3=1 

with strict inequality unless k = a. This completes the proof. □ 
Hence, if g(r, •) = for \r\ > Mo, then 

lim V" G(x)p n! N,M(x) = En !N;Mo (a;G). 

M— >oo ^ — ' 

From this it follows that the correlation kernel K(r, x; s, y) on the cylinder Z x Zyv 
is given by K(a; r, x\ s, y). We obtain the following proposition. 

Proposition 2.6. The correlation function for n = 2v + 1 non-intersecting walks 
on the infinite cylinder Z x Zjv as defined above is given by 



2nij /N \s—r ^27rij {x—y) /N 



(2.38) K{r,x;s,y) = ^J2(q+pe 2 

3=—" 

N-v-1 

s L (q+pe 2m il N ) s - r e 2m 3 {x ~ v)/N 



N-u-l 

1 

w 

N 

J=-v 

where LU r ^ s = 1 if r < s, uj r , s = if r > s. 
The induced measure on Z^r is given by 



■idet(i<r(0,^;0,^)) 1 < M)i/ < n =idet(i V e 2 ^'^-^)^ „ 

Til n I iv ^— ' 



1 J~J | e 27TM fl /Af _ e 27ria: l ,/JV|2 



l</i<^<n 

the equilibrium measure on Z^v (discrete CUE) , sec [^3| . 

We can take the limit n,N — > oo, n/N — ► p, < p < 1, and obtain a limiting 
dctcrminantal process on Z 2 . 

Proposition 2.7. TTie correlation function for the determinantal process on Z 2 
induced by non-intersecting random walks as defined above is given by 

(2.39) K{r,x;s,y)= (q + pe 27rze ) s ' r e 2me( - x - y) de 

J-p/2 

if r > s, and 

rl-p/2 

(2.40) K(r,x;s,y) = - {q+pe 2 ™ e ) s - r e 2 ™ e{x - y) d6 



p/2 



if r < s for < p < 1 . 
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This kernel is related to the B ± -kernels in f28|| . Compare also with 
3. Multi-layer discrete PNG 



We will discuss how the PNG model defined by (1.1), in the case when oj(x,t), 
(x, 4) e Z x N, satisfies uj(x, 4) = if 4 — x is even or if |x| > 4, can be embedded as 
the top curve in a multi- layer process given by a family of non-intersecting paths. 
We think of the u>(x,t):s as given numbers. The initial condition is h(x,0) = 0, 
x e Z. We extend h(x, 4) to all x e E by letting h(x, 4) — h([x], 4), which makes it 
right continuous at the jumps. Note that it follows immediately that h(x, 4) = if 
x < — 4 + 1 or x > t. 

For 4 — x odd we define the jumps, 

(3.1) r) + (x,t)=h(x,t)-h(x-l ) t) 

r~i~ (x, t) = h(x, 4) — h(x + 1, 4). 

We will see below that r) + ,rj~ > and we should think of r] + (x,t) as a positive 
jump at x at time 4, and 7/~ (x, t) as a the size of a negative jump at x + 1 at time 
t. Define 

(3.2) Tto(x,t) =min(r] + (x + l,<- 1), r)~(x- l,t- 1)) 
if t — x is odd and Tu(x, t) = if t — x is even. 

Claim 3.1. TTie jumps rj+ and rf~ satisfy the following evolution equation 

(3.3) r) + (x + l,t+ 1) = max(r/ + (a; + 2,t) - ?T (x, t), 0) + w(z + 1,< - 1) 
n + (x + l,t + 1) = max(?7"(a;,<) - n + (x + 2,t),0)+u(x + l,t- 1) 

/or i — x odd. Furthermore rj + (x, t) and T)~ [x, t) are > 0. 

Proof. We proceed by induction on i. Assume that ?7 + (x, t), rf~ (x, t) > for all 
x such that i — x is odd. We will prove that then ( |3.3[ ) holds, and hence 77+ (x + 
l,t + l),?y _ (x + 1,4+1) > for all x such that t — x is odd. Obviously our 
induction assumption is true for t = 0. Note that h[x + l,t) = h(x,t) — rj~(x,t), 
h(x + 2,t) = h(x,t) +n+(x + 2,t) - n~(x,t) and h(as- = h(a;,t) - n+(x,t). It 
follows from (1.1), our induction aasumption and w(x, i + 1) = 0, that 

ft(x + l,t + 1) = ft,(x,t) + max(0,?7 + (x + 2,4) - rj~(x,t)) + w(x + 1,4+1) 

/i(x,4+ 1) = /i(x,4) 

and the first half of fl3.3| ) follows. The proof of the second half is analogous. □ 

There is also an inverse recursion formula. 
Claim 3.2. If t — x is odd, then 

lu{x + 1,4+ 1) = min(r?"(x + l,4+ 1), ?7 + (x + l, 4 + 1)) 

(3.4) il + (x,t) = n+{x -1,4+1) -u(x- 1,4+1)+Tw(x -1,4+1) 
77- (x, 4) = 77+ (x + 1, 4 + 1) - u(x + 1, 4 + 1) + Tuj(x + 1,4+1) 



Proof. The first equation folows immediately from (|3.3|). F rom (3.2) and (^3) we 
see that the right hand side of the second equation in (|3.4| ) equals 

max(?7 + (x, 4) — n~ (x — 2, 4), 0) + mm(n + (x, 4), (x — 2, 4)), 

which equals ?7 + (x,4). The proof of the last equation is similar. □ 
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From this claim we immediately deduce the following 

Claim 3.3. If we know rj + (x + l,i + 1), rf~ (x + 1, t + 1) for all x such that t — x 
is odd, and Tu>(x, s) for s < t + 1 and all x, we can reconstruct u)(x, s), s < t + 1, 
a; £ Z, uniquely. 

Let a coordinate system be related to the (x,t) coordinate system via the 
transformation 



(3.5) 



(x,t) = (i- + j - 1), 



and define w(i,j) by (1.2). Then w(i,j) = if (i,j) r i? + , and this condition 
corresponds exactly to our assumptions on u>(x,t). Similarly to (1.2) we define 

(3.6) Tw(i,j)=Tuj(i-j,i + j-l). 

Claim 3.4. Assume that w(i,j) — if i or j is < s. Then, Tw(i,j) = if i or j 
is < s + 1 . 

Proof. It follows from the condition on w(i,j) that uj(x, t) = if (t + x + l)/2 < s 
or (i — x + l)/2 < s, whic h im plies, using (1.1), that h(x,t) = under the same 
conditions. It follows from (3.1), (3.2) and (3.6) that Tw(i,j) = if h(i — j + + 
j - 2) = or h(i - j - 1, i + j - 2) = 0. Now, h(i - j + 1, i + j - 2) = if i < s or 
j < s + l, and h(i — j — 1, i + j — 2) = if i < s + 1 or j < s. Hence Tw(i,j) = if 
i<s + lorj<s + l. □ 



It follows from claim 3.4 that T n w(i,j) — if i or j is < n, since w(i,j) = if 
i or j is < 0. Hence T n (ui(x, t) — if t < 2n — 1, since i+j — 1 <2n — 1 implies 
that i or j is < n. We formulate this as our next claim. 

Claim 3.5. Ift<2n-1, then T n uj(x, t) = 0. 

Let hi(x, t), i > 0, be the PNG process defined by ( pTl| ) with a>(x, t + 1) replaced 
by T l uj(x,t + 1), and with initial condition hi(x,o) = —i. We let T u(x,t + 1) = 
wfa;, t+ 1), so ho(x, t) = h(x,t) is our original growth process. It follows from claim 
|3.5| that at time t = 2n— 1 only ho, . . . , h n -\ can be non-trivial, i.e. hi(x, 2n— 1) = 
— i for all a; if i > n. Combining claim 3.3 and claim 3.5 we get 



Claim 3.6. Given hi(x, 2n— 1), x G Z, i = 0, . 

{w(x,i) ; f < 2n- l,x G Z}. 



. ,n— 1, we can uniquely reconstruct 



We can think of /ij at time 2n — 1 as a directed path from (— 2n + 1, — i) to 
(2rt — 1,—i) which has up-steps r] + (2m,2n — 1) at even ^-coordinates, x = 2m 
and down-steps r^(2ra,2n — 1) at odd x-c oord inates, x — 2m + 1, \m\ < n, and 
horizontal steps in between. According to 3.6 there is a bijection between these 
paths h , ... , h n -i and the set {ui(x, t) ; t < 2n — 1 , x € Z}. We set hi(x, t) = 
hi([x], t) for The paths obtained are nonintersecting: 

Claim 3.7. If t — x is odd, then 

(3.7) h l+1 (x,t) < hi(x-0,t) 
and if t — x is even 

(3.8) hi + i(x - 0,t) < hi(x,t), 
so that corners will not meet. 
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Proof. This is proved by induction on t. It is clearly true for t — 0. If it is true 
at time t it is still true after forming the maximum in (1.1) (deterministic step). 
(Note that hi and /ij+i have up-steps/down-steps at the same positions.) From the 
definition (3.2) it is still true after adding T l+1 uj{x,t) to the lower curve. □ 



In order to understand how a geometric distribution (1.3) on the w(i,j) is trans- 
ported to a measure on the non-imtersecting paths, we will assign weights to the 
jumps. Let and bj be given variables. The jumps are assigned weights as fol- 
lows: ri + (x,t) has weight a? i = (t + x+ l)/2 and T]~(x,t) has weight b v - ^ x,t \ 

j = (t — x + l)/2. Also, to T k u>(x, t) we assign the weight (aj, bj) T with the 

same correspondence between and (x,t), k > 0. The proof of the next claim 
is a straightforward computation using the defintions of the quantities involved. 

Claim 3.8. The product of the weights of rj + (x — 1, t + 1), ri~ (x — 1, t + 1) and 
Tlo(x — l,t + 1) equals the product of the weights of r)~ [x — 2,t), i] + (x,t) and 
lj{x - l,i + 1). 

Using this claim we can show that the measure is transported in the way we 
want. 

Claim 3.9. The product of all the weights of all the jumps in the multi-layer PNG, 
ho, . . . , h n -i, at time t = 2n — 1 equals, 

II M,r (lj) - 

i+j<2n 

Proof. Using claim [T^ repeatedly we see that 

II M;r (iJ) = II {a {t+x+l)/2 b(t-x + \)/2y^ 

i+j<2n xe1,t<2n-l 

= n ^ mM - x) K- { i mM - x) n (a it+x+1)/2b (t-x + i)/2f^ 

\m\<n x£Z,t<2n-l 



Repeated use of this identity toghether with claim 3.5 proves the claim. □ 



It is now easy to see that (1.5) holds. 
Proposition 3.10. Set G(i,j) = h(i — j,i + j — 1). Then 
(3.9) G(i,j)=m a x((g(i-lJ),G(i,j-l))+w(i,j) 
for i,j > 1. 
Proof. We have that 
h(i-j,i + j - 1) 

= max(h(i 2),h(i -j,i + j- 2), h[i - j + l,i+j - 2) + w(i,j) 

= max(G(i - 1, j), h(i - j,i + j - 2),G(i,j - l) + w(i,j). 

Since h(i - j - 1, i + j - 2) - h(i - j, i + j - 2) = r]~(i - j - 1, i + j - 2) > 0, this 
last expression equals the right hand side of fl3.9|) and we are done. □ 



If rj^,t]^ are the jumps for h r it follows from the assignments of weights that 
vfc(2m,2n — 1) = u has weight a^ + „ and r\~(2m,2n — 1) = u has weight b™_ m , 
\m\ < n, < r < n. If we think of the weights as labels transported from the 
w(i,j):s we see that if w(i,j) = for i > n or j > n, we have no labels a% with i > n 
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or bj with j > n and hence r]~(2m, 2n — 1) = if m < and 77+ (2m, 2n — 1) = if 
to > 0. Hence all plus-steps thake place to the left of the origin and all minus-steps 
to the right of the origin. This is the case discussed in EG] . From this consideration 



and (1.4) we obtain. 

Proposition 3.11. If\K\ < N, then 

(3.10) G(N + K,N-K) = h(2K,2N- 1). 
Also, ifw(i,j) = for \i\ > A or j > N, then for < K < N, 

(3.11) G(iV — K,N) = h(—2K, 2N — 1) 
and 

(3.12) G(N,N-K) = h{2K,2N-l). 

The discussion of the multi-layer extension of the PNG-growth model discussed 
above is closely related to the Viennot/matrix-ball construction, ||(J, [^i|, |^5| , of 
the Robinson-Schensted-Knuth (RSK) correspondence. Let us briefly discuss the 
relation. We can think of ( |3.2| ) and ( pT^ ) geometrically as follows. From (x, t) to 
{x— 1, t+ 1) we draw a line with multiplicity n + (x, t) and from (x,t) to (x+ 1, t+ 1) 
we draw a line with multiplicity r)~(x,t). A line with multiplicity zero means no 
line. At (x,t) a line with multiplicity rj + {x + l,i — 1) and a line with multiplicity 
r]~ (x — 1, t— 1) meet and we have a collision/annihilation of size Ttv(x, t) as given by 
(fj). If ry+(s + l,t-l) > r/-(x-l,t-l),thenr/+(a; + l,<-l)-?7-(2:-l,i-l)plus- 
lines survive and we add oj{x, t) new lines. Similarly in the other case. This explains 
(pT^). Assume that w(i,j) = if i or j is > A. If (w(i, j))i<i.j<N is a permutation 
matrix this gives exactly the "shadow lines" of the Viennot construction. We obtain 
a mapping to a pair of semi-standard Young tableaux P and Q of shape A. The 
number of to:s in the first row of P equals rj~ (—(A — to), A + to— l),m = l,...,A 
and the number of to:s in the first column of Q equals r/ + (—(A — to), A + to — 1), 
to = 1, . . . , A. Similarly, the same procedure starting with Tu instead gives the 
second rows and so on. Using this line of argument we obtain 

Proposition 3.12. Let (w(i, j))i<ij<N be given and set w(i,j) = if i or j is > 
A. The RSK- correspondence maps a submatrix (w(i, j))i<i<M,i<j<N > M < N to a 
pair of semi-standard Young tableaux of shape A(M, A) = (Ai(M, A), Aa(-M, A), . . . ). 
(Similarly, we can consider (w(i, j))i<i<N,i<j<M ■) Consider the family of height 
curves hi, < i < N, obtained from the multi-layer PNG process using (w(i,j)). 
Then, forO<K<N,l<j<N, 

(3.13) Xj(N-K,N) = ^_i(-2A,2A- 1) + j - 1 
and 

(3.14) Xj(N, A -K) = hj-x{2K, 2A - 1) + j - 1. 

If we add vertical line segments to the graphs, x — > hi(x, 2 A — 1), < i < A, we 
obtain A non-intersecting paths with hi(-(2N- 1), 2A- 1) = hi(2N- 1, 2A- 1) = 
— i. Recall that hi(x, 2A— 1) = 1 — i for i > A so that at most A paths are "active" . 
The paths are described by particle configurations. Let 

C2N-x{x) = (h (x, 2A - 1), . . . , h N -i{x, 2N - 1)) 

and 

CW-i = (C 2N -i(-M + 1), . . . , C 2N -i{M - 1)), 
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where M = 2N - 1. Note that CW-i(-M) = CW-i(M) = (0, -1, . . . , -N + 1). 
Set 



(3.15) 
(3.16) 



<hi-i,2j[x,y) = \ n J+ 

it y < a;, 







if y > x, 

v (1 - b N - :j )b x N y. J \iy<x 

for |j| < iV with the convention that 0° = 1. It follows from the Lindstrom- 
Gessel-Viennot method or from the Karlin-McGregor theorem that the weight of 
the non-intersecting path configuration corresponding to C 2 n—i — x, with weights 
assigned to jumps as above, equals 



The way the weights are related to the "weights" of the geometric random variables 
as described above shows that 

M-i 

P[C 2N _ 1 =x] = - [] det(0 r , r+1 (<,^ +1 ))^ =1 



,M 



r=-M 



with Z n ^M g iven by ( 1.1 8| ) and n = N, M — 2N — 1. Hence we obtain a measure 
of the form ( 1.17 ) . We note that 



nSi(i-«i) n (i-&i) f 



(3.17) 



Z n ,M — 



11, r v/0 

We summarize what we have found in the next proposition. 

Proposition 3.13. Let hi, i > 0, be the multi-layer PNG process obtained from 
geometric random variables with parameters ciibj as defined above and let <fir,r+i be 
defined by (SAL) and (3.16). Then, 

(3.18) F[h k -!(r,2N- 1) =4,1 < i < n, \r\ < M] 

M-l 

— J] det(^ r+1 (<,^ +1 ))^ =1 , 



Z .I 



A I 



r=-M 



where Z n ,M * s given by (3.17), x i — xf 1 = 1 — i, x\ > x 2 > • •• > 2;^ /or eac/i r, 
n = N and M = 2N — 1. 

The fact that the probability measure has this form makes it possible to compute 
the correlation functions. Set 



(3.19) 

and 
(3.20) 



/„_!(*) = (1 - a j+N ) °?+NZ m = 



1-a 



j+N 



f 2j (z) = (l-b N - j )J2b%- j z m 



1 - a j+N z 

1 - b N _j 
1 - bN-j/z 



so that (1.25) holds. The interpretation of the correlation functions given by (1.23) 
in this case is that they give the probability of finding particles at the specified 
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positions. We can take n > N in ( 3.18 ), where n is the number of PNG height 
curves. All height curves hi with i > N have to be trivial, i.e. hi = —i if i > N. It 
follows that the probability of a certain configuration is independent of n for n > N . 
Thus, we can take the kernel K n ' M , (1.22), with an arbitrary n, n > N arbitrary 
and obtain the same value. In particular we can let n — > oo and use proposition 
1.8. It is clear that all the conditions of this theorem are satisfied when f r (z) is 
given by ( | 19[) and fl3.20| ). Let r = 2u, s = 2v both be even, |u|, \v\ < N. The 
expression (1.3C) becomes 



(3.21) 



G(z, w) = 



Xlj=u (l- 



EC 



/ 1 — q]V4 
-N+l \ l-a N+j 



+3 
I W 



llj=-iV+l [l-a N+j /z J llj=v [l- 



— bj^—jW 



We summarize our results for the correlation functions in a theorem. 

Theorem 3.14. Let the multi-layer PNG process be defined using geometric ran- 
dom variables w(i,j) with parameter aibj, < a,i,bj < 1, and let G(z,w) be given 
by ( POO) ;. Set 

1 



(3.22) 



K N (2u,x;2v,y) 



(2niy 



dz 
z 



dw w y 



-G{z,w), 



where "f r is the circle with radius r and center at the origin, 1 — e <r\ <r2<l + e 
with 1 + e < min(l/&j-), 1 — e > max(aj) and \u\, \v\ < N, x, y € Z. Furthermore, 
let 

huM^V) = ^ f e l ^ e G(e w ,e l9 )de, 



(3.23) 

for u < v and 4>2u,2v{x, y) = for u > v. Set 



(3.24) 
Then, 
(3.25) 



K N (2u, x; 2v, y) = K N (2u, x; 2v, y) - (f>2 U flv(x, y). 



V[{2u,xf ) G {(2t,hi(2t,2N - 1)) ; \t\ < N , < i < N}, \u\ < N, 1 < j < k u ] 
= det(K N (2u,xf u ; 2v,x 2 j v ))\ u ^ v \ <N:1 < l <k u ,i<j<kv 
for any x 2u G Z and any k u € {0, . . . , N}. 

Consider the finite-dimensional distribution of ho(x,t) — h{x,t), the top curve, 
(3.26) W nM [h (2si,2N - 1) <£ t ,l<i< to], 

where n is the number of paths, M = 2N - 1, n > N, |s, | < N and li > -N. This 
can also be written 

Pn,Af[no particles in {2si} x (li, oo), 1 < i < to]. 

This probability is independent of n > N and hence we can let n — > oo. Let 
g(2si,x) = — Xti,oo) ( x )i 1 ^ i < to, and g(r,x) = if r is not equal to one of the 



2sj. Hence, by proposition 24 

(3.27) F N M[h (2s t ,2N - 1) < ^,1 < * < to] = det(7 + gK N ) L 2 {Ku) . 

This formula can be used to study the convergence in distribution of the rescaled 
height curve. 
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4. ASYMPTOTICS 



We will consider the asymptotics of the kernel ( 3.24 ) in the case a« = b% — a for 
all i > 1, so that w(i,j) are geometric random variables with parameter q = a 2 . 
The function G(z,w) in (3.21) then becomes 

(4 1) G(z w) = (l- q)*-«> (1 ~ a/z)N+U (1 ~ aw)N ~ V 

Write 

1 (z - a) N+u 



z x+N+u (l- az )N-u> 



so that, by ( |3.22[ ), 

(4.2) K N (2u,x;2v,y = V ; / — / 

(2tti)^ J lr2 z J lri w z-w w 

where a < r\ < ri < 1/a. 

Set fi = m/N, fx' = m'/N, (3 = u/N, (3' = -v/N and 

UA Z ) = j^\ogF u , m _ N (z) 

= (l+(3) log(z -a) -(1-0) log(l - az) -(n + 0) log z. 

Then, f'(z) = P(z)/Q(z), where Q(z) = z(z — a)(l — az) and 

n/ a / a\r 2 1 - a 2 - /i(l + a 2 ) M + /9i 

a(/x — p) jtt — p 

We will write 

(4.3) p = pOhfi = ; i = a ^V = —p- 



The critical points of / are z^r — p ± \/p 2 — q an d we obtain a double critical point 
if p 2 = q which gives 



(4.4) M = Mc (/3) 



1 + a 2 //l + a 2 \ 2 4a 2 /3 2 



1-a 2 \l\l-a 2 J (1 - a 2 ) : 

and 



(4-5) p c = p{/j, c , /3) = 



1 + a 2 + 2a v / l~/3 2 - (3(1 - a 2 ) 
Set 

, . , cW^l + a) 1/3 „ 1 - a , 

(4.6) d= >— , d' = -—d. 

1 — a 1 + a 

It will be convenient to write 

(4.7) u= L tN V3 v = }_ T < N 2/3 j 

a' a' 

since N 2 / 3 is the right scale for u and v if we want a non-trivial limit. The correct 
way of writing x and y will turn out to be 

(4.8) x = N(n c {l3)-l)+ZdN 1 / 3 , y = N(^')-l)+^'dN 1 ^. 
We will assume that |r|, |r'|, |£|, are < log AT. 
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The paths of integration can be deformed into 
(4.9) T:R3t'^ z(t') = p e (fi) + - JjL^ = p it, 



(4.io) r':R3/. W (s') = Pc (p>) + -JL m --£ m ±y- it )-\ 

where 77, 77' > will be appropriately chosen; we will require that 
In that case we have, by Cauchy's theorem, 

4.12 K N (2u,x;2v,y) = - K > /— / F ux (z)F- vy (—). 

{2iri) z J r z J r i w z — w w 

We first estimate this integral and then we will compute its asymptotics using a 
saddle-point argument. Choose fi so that p(/j>, /?) = p = p c (P) -fry/cL/V 1 / 3 as in (4.9), 



and let q — q(p, (3) be the value we get with this [i. This is possible by formula 
Ne can -\ 

P'n.x (^) 



(4.5) with (i > [i c . We can write 

1 1 (z - a) N+u 



z x+N-nN Z {^i+I3)N (I _ az ~)N-u ' 

and then let z = p — it, t e R, and take the absolute value to get 
(4.13) \F U , X (P-U)\ 2 = {p2+t2 l )x+N ^ N e 2Nh ^\ 

where 

2h(t) = (1 + P) log A - (1 - 0) log 5 - (fx + 0) log C, 

with 

A = (p - a) 2 + t 2 = 1 - 2ap + a 2 + — ^- + p 2 - q + t 2 

fi-P 

B = (1 - apf + a 2 t 2 = 1 - 2ap + a 2 + ^JL + a 2 (p 2 -q + t 2 ) 

C = p 2 + t 2 = l + ^-+p 2 -q + t 2 . 
fi-P 

Note that 

(A* - P)A = (1 - a 2 )p + 1 - a 2 + (/1 - P){p 2 - q + t 2 ) 
(p - P)B = 1 - a 2 - (1 - a 2 )p + a 2 (fi -fj)(p 2 -q + t 2 ) 

(p - P)C ={ t l-p){p 2 -q + t 2 )+n + P. 
A computation now gives 

(4-14) h\t) J^l±fl {(l - a 2 ) 2 - (1 - + P) + (1 + a^p + 2a 2 p 2 

-a 2 (n-P) 2 (p 2 -q + t 2 )}. 
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Another computation shows that p 2 — q > 0. To leading order we have p 2 — q « 
2r)/dN 1 / 3 . Recall that t = t'/dN 1 / 3 . When t is large we have h'(t) » -(/z - /3)/t, 
and we also have the estimate 

(4-15) h'(t) < (Af _ /3 t )Aj3C (p 2 -g){(l-a 2 ) 2 -(l- Q 4 )( / i+/?) + (l+a 4 )M/3+2a 2 /3 2 }. 

for i > 0. Consider the case < t' < N 1 ; the case — N 7 < t' < is analogous by 
symmetry. Here < 7 < 1/3. Using ( |4.15| ) we see that - /i(0) « -2r)t' 2 /N 
and we can show that 

(4.16) 

for < iV 7 , < 7 < 1/3, and AT sufficiently large. Define h*(t,p) by 



Nh'(t) < -^rjt' 2 + Nh{0) 



a 2Nh,(t,p) _ 



1 



[(p-a) 2 +t 2 ]^ N 



(p 2 + t 2 )^+P) N [(1 - pa) 2 + Q2 t 2](l-/3)JV ■ 

A computation, compare (4.23) below, gives 

e Mo,Pc) ^(i_ a )2« e ^ 3 

and K(0,p) = h*(0,p c )+d 3 (p-p c ) 3 /3+- ■ ■ = h*(0,p c )+r} 3 /3N+. . . . Consequently, 

1 (1 - a) 2u (T 3 



D JVh(0) 



P 



Combining this with (4.16) gives 



Write 



\F u , x (p-it)\ < 



1 



C 



(p 2 + i2)( :C +JV(l-p))/2 p N(p-^) 



(1-q) 2 " (r 3 +T> 3 )/3 _ 3 ^2 /2 



(p 2 + £ 2 )(z+A r (l-M))/2 

Further computation shows that 



{p-it) x + N 0--Mc) 



p 2 +t 2\N(n-fi c )/2 



(p-uy+xi 1 -^) 



(l + * 2 /i> 2 ) JV( "~" c)/2 ^e" 4 ' 2 
Collecting the estimates we find 



ITs)! <C(l-a)" 



2u e ^(T 3 W)-ir-i V -^t' 



(4.17) 1^^) + ^ ^ 

for \t'\ < AT?, < r) < A"?, < 7 < 1/3. 

Using (4.14), (4.15) and the other estimates above we see that the contribution to 
the integral from \t'\> N~< and/or \s'\ > N" 1 is < Cexp(— cN 2 " 1 ) for some constant 
c > 0. Hence, using the parametrization ( |4.9| ) in (4.12) we can restrict to < N" 1 , 
\s'\ < N 1 . We can use ( 4 . 1 7| ) if we want an estimate of the integral. To get the 
asymptotics we make a local saddle-point argument. 

To leading order we have p c {/3) = 1 + r/dN 1 / 3 , p c (@') = 1 - T'/dN 1 / 3 and hence 
the condition (4.11) requires 

(4.18) r-r' + rj + r?' > 0. 
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We will use the parametrizations (4.E) and consider the integral 
(4.19) 



(1 - q) 2 ("-") 
(2m) 2 



dt 



ds 



z'(t) w'{t) z(t) w ( s )y+N(i-MP')) 
,\ a ./M<arr ■""*(*) w(t) z(t)-w(s) z(t)*+W-i"W) 



Now, 



(4.20) Nf^pipM + j^iv - H)) 

where the remainder term rjv(i) can be neglected for |t| < JV 7 . Also, z'(i) 
-i/dN 1 ' 3 , w'(s) = iw(s) 2 /dN 1 / 3 and 

(4.21) »« diVV3 



Furthermore 
(4.22) 



z(t) — w(s) t' — t + i(t + irj + s + irj') 

w K b ) _ £'T'-£T+i£{t+in)i£'(s+iTi') 



We also need to compute 

m (pM-a) N +u 



1 



(l-ap c ((3)) N - u P ?^ {P} - U ' 



Using the formulas (4.4) and (4.5) above a rather long computation, which we omit, 
shows that 

(4 23) v r ' ; ' 1 } ' " 



,NUcW),e(Pc(P)) ~ (1 - a) 2u e r3/3 . 



Inserting (|4.20| ) - fl4.23|) into fl4.19| ) we see that, provided fl4.18| ) holds, 
(4.24) 

lim dNV 3 K N (2 — -N 2 l 3 r, —N + (£ - r 2 )^ 1 / 3 ; 
Af-too 1 — a a 1 — a 



1 + a 1 
1-ad 



^2/3^ 2a - ■ <+> > 2 - 



1 — a 



N + (£' — t )dN 1 / 3 ) 



47T 263 



(t 3 -t' 3 )+C't'-Ct 



lvaz—rj Jlmw—r/ 



+ i(z + w) 



-dzdw. 



Here we have used /ti(/3 c ) = ~ 1 i > ct 2 /3 2 + We also want to compute the 

corresponding limit of (3.23) with G(w,w) given by fl4.1|), i.e. we consider, u < v, 



(4.25) 



2ir /__ 



If we set g(8) = log(l + a 2 — 2a cos 8), then 



g'(0) 



2a sin 9 



1 + a 2 - 2acos9 : 
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and we see that g(9) has a quadratic minimum at 9 = 0. Hence, we can immediately 
both compute the asympotics of and estimate the integral in ( 4.25 ) when x — 2a(l — 
a^N + (£ - T 2 )dN^ 3 , y = 2a(l - a)" 1 ^ + (f - r'^dN 1 / 3 , u = ^d~ l N 2 / 3 T 
and v = ^d^N 2 ' 3 ^. We obtain 



(4.26) 



lim dN'^fouMx^y) = ^~ I e l(i ' 



-i+T 2 -T' 2 )t-(r'-r)t 2 dt 



\J^(j' - T) 

We want to identify the right hand side of (4.24) combined with ( 4.26| ) with the 
extended Airy kernel. This can be done using the proposition 2.3. Combining 
this double integral formula for the extended Airy kernel with (1.24) and ( 4.26| ) we 
obtain the following result. 



Proposition 4.1. Let d — (1 — a) 1 cr L / 3 (l + a) 1 / 3 and let K N be given by ( 3.24 )■ 
Then 

lim dN^K^l^d^N^T, —A +(£- T 2 )dN 2 / 3 ; 
1 — a 1 — a 



N— >oo 



1 



— n + (e - T> 2 )dN 2 ' 3 ) 



(4.27) 



- Jt 3 -t' 3 )/3+Z't'-Zt 



uniformly for £, r, r' in a compact set. 



We can now combine the formula 3.27, theorem 3.14, proposition 4.1 and some 
estimates of Kn, which can be obtained from the asymptotic analysis above, to 
prove the following theorem on convergence in distribution to the Airy process. A 
complete proof requires a control of the convergence of the Fredholm expansions 
but we will not present the details. The individual determinants in the Fredholm 
expansion can be estimated using the Hadamard inequality. Compare with theorem 
1.2 and lemma 3.1 in [BTstl. 



5. A FUNCTIONAL LIMIT THEOREM 

5.1. A moment estimate. Consider the PNG height functions hk(x,2N — 1) 
defined in sect. 3. Set 

t ■ = — j — 

h cN 2/3 

where c = (1 + a)(l — a)~ 1 d~ 1 , j G Z. The normalized height functions are 

1 2a 



with d as in (p~q), k G N. For a given function / : K 

1 , 2a 



we write 



Mx) = /(- 



■TV)). 



'dN 1 / 3 ^ I- a 
Assume that there is a K such that f(x) = for x < K. Define 



(5.1) 



H N {f, tj) = J2 MfoW, 2N -!))=£ f{HN,k(tj))- 



fe=0 
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Lemma 5.1. Assume that f is a C°° function and that there are constants K\ and 
K 2 such that f{x) = if x < K\, and f{x) equals a constant if x > K 2 . There is 
a constant C(f,a) so that 

(5.2) n(H N (f,t u ) - H N (f,t v )) 4 } < C(f,a)e~^ 3 \t u - t v \ 2 , 
for \t u - t v \ < 1 and \t u \, \t v \ < \ogN. 

Proof. The proof is rather long and complicated. We will outline the main ideas 
and steps in the argument without giving full details. The left hand side of (5.2) 
can be written 

00 4 

(5.3) J2 nY[(fN(h kr (2u,2N^l))- f N (h kr (2v,2N-l)))}. 

ki,k2,k3,k^ — l r — 1 

We can rewrite this using formula ( 3.25| ) in theorem 3.14. Let us write the kernel 
Kn(2u, x;2v,y) in (3.24) as K uv (x,y). We will use the following notation: 

(5.4) K\ --I; ^ ^ ^ = det(A-(r ijiBi ; aj>W ))& =ll 
and we will also write 

(5.5) K(rx,xi r 2l x 2 ... r m , x m ) = det(K (n, x^, r j} Xj))^ j=1 . 
Furthermore, we will write 

(5.6) D ) = K(2ui,xi 2u 2 , x 2 ■ ■ ■ 2u m , Xm). 

Set h N (xi,x 2 ,x 3 ) = -6[fN(xi) 2 f N (x 2 )f N (x 3 ) + fN(xi)f N (x 2 ) 2 f N (x 3 ) 
~fN{xi)fN(x 2 )fN(x 3 )], which is synmmctric under permutation of x\ and x 2 . 
Then, the sum in (5.3) can be written 

(5.7) 

^ fN(xi)fN(x 2 )fN(x3)fN(x4)[D uuuu (x 1 ,X 2 ,X 3 ,X4) ~ 4D uuuv (xi, X 2 , X 3 , X 4 ) 

zez 4 

+ GDuuw{xi,x 2 ,x 3 ,x i ) - 4D uvvv (xi,x 2 ,x 3l x i ) + D vvvv (xi, x 2 , x 3 , X4)] 
+ { 6 ^ JV (^i ) 2 ^ w (^2 )/tv (x 3 ) [D uuu (xi , x 2 , x 3 ) + D vvv (xi , x 2 , x 3 )] 

+ h N (xi 1 x 2 ,x 3 )D uuv (xi,x 2 ,x 3 ) + h N (x 3l x 2 ,xi)D uvv (xi,x 2l x 3 )} 

+ ^ 2 (fN(xi) 3 f N (x 2 ) + fN(xi)f N (x 3 ) 3 )[D uu (xi,x 2 ) - 2D uv (xi 1 x 2 ) + D vv (xt,x 2 )] 

+ 3/jv(xi) 2 /iv^2) 2 [^ nil 

(x!,x 2 ) + 2D uv (x 1 ,x 2 ) + 

DyV (X\ , X 2 ^j\ 

+ ^2 2f N (x 1 ) i [K uu {x ll xi) + K vv (X\ , X\ Jj 
= E 1 + E 2 + E 3 + E 4 + S5- 

We have K uv = K uv — <j) if u < v and K uv — K uv if u > v. Here we have written 
(friijV- Se t AK UV — K uv K uu , AK VU — K vu K uu ^K vv — Kuv K uu . 
We see from (4.26) that cf> acts like a kind of approximate ^-function. This will 
be important for the cancellation between different terms in ( |5.7| ). The argument 
goes as follows. We will take out all terms in (5.7) containing <f> and combine them 
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with other terms so that we get cancellation. We will then expand in AK UV , AK VU 
and AK VV . The terms linear in AK will cancel and what will remain will be terms 
containg (AK) 2 or higher powers. They will give a contribution proportional to 
\tu — t v \ 2 which is what we want. 

In the computations below we use symmetries and also relabelling of variables. 
Expand in <f> and in the terms linear in <f> we expand in AK. Let D denote the 



same object as in ( p.Q ) but with K replaced by K. We find 
(5.8) 

Si = !N(xi)f N (x 2 )f N (x 3 )f N (xi)[b uuuu (xi,X 2 , X 3 , Xi) - &D uuuv (xi,X2,X3,X4) 

x 

+ 6D uuvv (x%, x 2 , x 3 , Xi) - AD 

+ Y^2^{ Xl ,x A )KZ fa X4 ) [AK vu (x 2 ,x 3 ) + AK uv (x 2 ,x 3 )- 
x \ Xl X2 J 

AK vv (x 2 ,X 3 )]f N (xi)f N (x 2 )f N (x 3 )f N (x4) 

-^^(xi-,X4)(i) ( X2,X 3 )kll ^ } N {x 1 )f N (x 2 )fN{x 3 )j N {x i ) 

^2 (terms with AK 2 ). 

x 

We will give a brief discussion of the AK 2 terms below. Also we will see then that 
terms containing 

(5.9) AK vu (x,y) + AK uv (x,y) - AK vv (x,y) 

will give a contribution proportional to \t u — t v \ 2 . If we expand the _D-part of ( |5.8D 
in AK we will see that the terms linear in AK cancel out. Since obviously the 
0:th order term equals zero we are left with A_fT 2 -terms. The term containg two 
</>-factors will be combined with other terms below. 
We expand £2 similarly. The part linear in AK is 

(5.10) -Y / WN(x 1 ) 2 f N (x 2 )f N (x 3 )KZ fa ^[AK^Xs) 

+ AK uv (x 2 ,x 3 ) - AK vv (x 2 ,x 3 )}. 



Actually this sum can be combined with the corresponding term in (5.8) to get 
some cancellation, see the 0-calculations below, but we can also use the fact that 
( pi} ) has the right order. We also get Aif 2 -terms and a term linear in cf>, 

(5.11) 

^2l2<j)(x 1 ,x 3 )[f N (x 1 ) 2 f N (x 2 )f N (x 3 ) ~ fN{xi)fN(x 2 )f N {x 3 ) 2 ] 

X 

X \ K ™ [ Xl x 2 ) K ™ [ Xl x 2 

+ Y J V<}>(x l ,x 3 )f N (x l )f N (x 2 ) 2 f N {x 3 )i y Kz{^ i X ^)+Kl(^ x \ 



b- 
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Consider next S3. We get a term linear in (j), 

(5.12) -J^zi.ara (x 2 , Xl)[fN(xi)fN(x2) 3 + fN(xi) 3 fN(x2)] 

x 

a term linear in AK, 

(5.13) ^^(uJ/at^) 3 + f N (xi) 3 f N (x 2 )])[AK vu 

x 

AK uv (x u x 2 ) - AK vv (x u x 2 )]K uu (x 2l x 1 ) 1 

and Aif 2 -terms. In ( 5.13|) we again have the expression ( |5.9| ). 

The leading term in £4 is 
(5.14) 

Y / 3fN(xx) 2 f N (x 2 f 



K v , 



X\ x 2 



+ 2K™ Xl X2 

\Xl x 2 



K vv 



Xl X 2 
Xl X 2 



K Xl X 2 

and we also have a term linear in 

(5.15) ^6fN(xi) 2 fN(x 2 ) 2 (/)(xi,x 2 )K vu (x 2 ,xi). 

X 

Finally we have £5 which is 

(5.16) Y, VN{xif[k uu {xi,x x ) + K{xi,xi)]. 



When calculating the cancellations involving the </>-terms we will combine the double 
<^)-tcrm in (5.S) with £& in (5.11) and (5.14). Also we will combine (5.12), (5.15) and 
( 5.16] ). We will discuss this second case first in some detail and then the first case 
more briefly. The term S Q is similar and finally we will indicate what is involved 
in estimating (5.9) and the Aif 2 -terms. 
We want to estimate 

(5.17) Hx,y)K vu (y,x)[6f N (xi) 2 f N (x 2 ) 2 ~Af N ( Xl )f N (x 2 ) 3 

- 4fN(xi) 3 f N (x 2 )] + J2( R ^(y> y) + K™(v> y))fN(y) 4 - 

Here we have made a symmetrization in x and y by setting 
K vu (x,y) = -[K vu (x,y) + K vu (y,x)]; 

note that (j)(x,y) is symmetric in x and y. Next, we will introduce some notation 
and some formulas that will be used. Set 

1 

z 

and 



(5.18) 



(I-^ + -- 2 ) 



1 — a/z 



1 — az 



(5.19) G* ab (z,w) = 
so that 

(5.20) K ab (x,y) = 



N 



1 — aw 

1 — a/w 



N 



(l + g(z)r(l + g( W )f 



1 



w(z — w) 



(2ni) 



dz / dwG* ab (z,w) — , 
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1 



and 
(5.21) 

Note that 

(5.22) G* ab (z, w)(l + g(z)f{\ + g(w)y d = G* a+cMd (z, w). 

Fix e > and let f e {x) = f(x)e~ ex . Then / e is in L 1 (IR) and we have 



(5.23) 



f N (xY 



1 



(27T) 



F^{X)e i ^ {x ^ cN ^d m X, 



where F£(A) = / 6 (Ai) . . . M)mh c = 2a ^ + a)" 1 and £ m (A) = (X ± + ■ ■ ■ + A, 
ime )/dN 1 / 3 with d given by ( |4.8| ). Integration by parts gives 

(5.24) (1 + g(z)) u - v -l = (u- v)g(z) + -^^Ri(z) + WtF^W, 



d 4 7V 4 / 3 



where 



and 



Rx{z) = d i N 4 ' 3 g{zf 



1 - t 



o a + tg(z)) 



-dt 



R 2 (z)=d 4 N 4 / 3 (log(l+g(z))) 2 / (l-t)(l+g(z)) t ^dt 



Let A/(x) = f(x + 1) — /(x) be the usual finite difference operator. We have the 
following formula 

(5.25) 

v) + % ) M*) m = My) m 5 [(1 + ff(z)) "^ + (1 + 5H) n ^l 

a(it — i?) 



{l-af 



.,y+i 



T A/j?(y) 



— A/J?(y-1) + — A/ff(y) 



+^rA/^( y -i) 

J u — w 
t d 4 iV 4 / 3 

(u — v) 2 
+ d 4 N 4 / 3 



(27T) 



d m A J F 1 4(A)e l ^( A )^- cAr ) 
^(ze-^-W) - i?i(z)i?i(we l?m(A) ) - EiH 
i? 2 (ze-^< A )) - i? 2 (z) + i? 2 ( W e 4 «-"( A >) - R 2 (w) 



for | to | = exp(— me/dN 1 / 3 ) — r%, \z\ = r 2 = 1/ri. To prove this, introduce the 
formula ( 5.21 ) for and the formula ( 5.23 ) for fff into the left hand side of ( |5.25 ) 
and use 



1 

dNV* 



(Ai + --- + A m ) 
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where So is the Dirac (5-function, to carry out the ^-summation. This gives 
(5.26) 



1 



(27T) 



< r n XF'(A)e«"W&- cJV > — 



x [(1 + g(ze- tUW )) u - v + (1 + .g( W e^" (A) )) u - 1 '] 

= My) m ^[(i + g{z)) u - v + (i + jWH 

+ (i + g (we^W)r~ v - (i + g(w)r- v }. 



In the last expressi on w e use ( 5.24 ) and the explicit form ( 5.18Q of g to obta in th e 
right hand side of ( [5.25 ). We will call the first part of the right hand side of ( 5.25 ), 
/Af(?/) m 7^[(l + d{ z )) u ~ v + (1 + 9( w )) u ~ v ]i the contraction term, which is the main 
contribution. The second part is called the finite difference term. 

We can now insert the integral formula ( 5.20| ) into (5.17) and use ( |5.25 ). The 
contraction term from the first sum in (5.17) will then exactly cancel the second 
sum. Here we use (5.22). What remains is 



(5.27) 



u — v 

d 4 A 4 / 3 



So 



u-v (u - v) 2 

1A AT/1 ii b 2, 



d 4 A 4 / 3 



^4^4/3 



where So is the part coming from the finite differences, S± is the part coming from 
terms involving R\ and S% from the terms involving i?2. After some computation 
we find 



(5.28) 



Also, 
(5.29) 



S Q = - 



(1 - a) 2 (2tti) 2 



7r; 



<l: / dwG* vu {z,w){ — hro) 



*E^(^ 1/3 (M»+i)-Mi/))) 4 - 



3 (A)(v-cJV)^_ 



1 



(2ttz) 2 



dz I dwG* vu (z,w)— {6[h(z,w;&W) 



-h{z, w; 0)] - A[h(z, w; 6(A)) - h(z, w; 0)] - i[h{z, w; 6(A)) - h(z, w; 0)]} , 

i = 1,2. In order to restrict the y-summation so that (y — cAQ/d/V 1 ' 3 ranges 
over a compact interval we make a summation by parts in ( |5.29 ). Recall that we 
assume that f(y) is a constant for large y. If we let a = eKp(i^(X))w / z and use 
a y = (1 — a)^ 1 (a y — a v+1 ) in the y-sum in ( 5.29| ) a summation by parts gives 



1 - exp(i^ 3 (X))w/z 



zv 
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Hence, for i = 1, 2, 



S - - Mi/ - 1))7233 /. d3AF 3 

Iff w v 

/ 27r ^ a y y d"wG vu {z,w)— {6[h(z,w;&(X)) 

- h(z, w; 0)] - 4[/i(z, w; £ X (A)) - w; 0)]- 
4[fe(*,u>;6(A)) <))]}■ 



^?3(A)( a -cN) 



1 - exp(j£ 3 (A))w/z 

The expressions Sj will be estimated using the types of estimates derived in sect. 
4. 

Write u = (1 + a)(l - a^cHiV 2 / 3 -/-, v = (1 + a)(l - a)" 1 ^ 1 ^ 2 / 3 ^ and 
y = 2a(l — a) _1 iV(^ — T 2 )dN 2 ^ 3 . To estimate (5.28) we can now use our results 
from section 4. We use 

ri it t + n — it 

^ ^=^ + ^-^ = 1 + ^7^ + --- 

as parametrizations of the integrals as before. Using the same estimates as in section 
4 we can restrict the integration to \t\, \s\ < with an error < Cexp(— cN 2 " 1 ) 
with some c > 0. Since v — u > 1, and hence t v — t u > c/N 2 ^ 3 , and fu rthe rmore 
| t tt | < log A T, we can incorporate the error term into the right hand side of (^^) . The 
integral in ( [5.28 ) can then be estimated using ( [4.17 ). Note that, by our assumptions 
on /, the number of y-terms ^ is < CN 2 / 3 and we get a compensating factor 
1/N 2 / 3 from the parametrizations; see also (4.21). The numbers 77,77' are chosen 
so that T) + 7/ > 2. Note that, since we assume r' — r < 1, the condition ( 4.18j ) is 
satisfied. We find 

(5.32) \S \ < C(f, C*)-^ 2 e (r 3 -r' 3 )/3+ S (r'-r)+(^+V 3 )/3-«( J7+ V) ) 

2/ 

where the y-summation is over ally G Z such that [y-cN) / dN 1 / 3 e [Jfi-l,Jf 2 +l]. 

Consider now Sj. Write z = z exp(-me/diV 1 / 3 ), A = (Ai + h X m )/dN 1 / 3 . 

Then, 

g^e-^W) = - " p- i) 2 l e ^ - 2i(z- l)sinA + 2(cosA- 1)]. 
(1 — a) z 

Thus, 

1 + 9 (ze-^ A >) = 1 + —^(1 - cos A) - - - \^—^e fx - 2i(z - 1) sin A]. 
(1 — ay (1 — a) A z 

The last term is small for large N and the second is > 1. Hence, 
(5 - 33) \l + tg(z(t)e-^W)\ ~ 2 

for \t\ < N 1 and N sufficiently large. Consequently, there is a constant c%(a) 
depending only on a such that 

(5.34) | 3 (z e -^^)|< Cl (a)(A 2 + |z-l| 2 ) 
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and 

(5.35) | log(l + tg(ze-*«>M))\ < ci(a)(A 2 + \z- 1| 2 ). 

We can also write 



(5.36) 



1 + tg(ze 



-«m(A)i 



(1 



2a 



r(l-cosA)) 



1 



{l-af 
{z-l) 2 e il /z-2i{z-l)sm\ 



(1-a) 2 2a(l-a)- 2 (l-cosA) 

By periodicity it is enough to consider |A| < tt. Estimating the cosine and sine 
functions we see that there are constants 02(0;) and 03(a) such that 

|1 + tg(ze~ l ^)\ > exp(c 2 (a)A 2 - c 3 (a)(\z - 1| 2 + \z- 1||A|)) 

and hence, 

|1 + i s (*e-*™W)|*< tt - > < exp(-c 2 ( a )A 2 + c 3 (a)(\z - 1| 2 + \z - 1||A|)). 

Estimating the quadratic polynomial in A we obtain 

|1 + tg(ze- iUW )\ tl - u - v) < exp(t(r' - T)c 4 (a)N 2/3 \z - 1| 2 ). 

Now, z(t) exp(— me/dN 1 / 3 ) = 1 + (r + rj — me — i^/dN 1 / 3 + . . . , and we obtain an 
estimate 

(5.37) |1 + tg(z{t)e~ l ^)\ t{ - u -^ < exp(c 5 (a)[(r + r? - me) 2 + t 2 }). 

A computation shows that 

1 



(5.38) 



< c 6 



|1 - e ^3(A)«i/z| 

if r — t' + rj + if > e > 0. Since t't < 1 and we take 77 + 77' > 2 we see that we can 
take e = 1/2 for example. Furthermore, since (y — cN) / dN 1 / 3 is bounded for the 
y:s that contribute to the sum, 

(5.39) | e ^ 3 (A)(y- c Jv)| ^ ct _ 

We can now again estimate as in sect. 4 and use (4.17). This results in an estimate 

(5.40) \Si\ < c s (t,a) (J(l + \*)\r(\)\d\\ e M a ^ + v- m ^ 



N 1 / 



1/3 



)/3+{(r'-r) + (^ 3 +r ) ' 3 )/3-5(r,+» ) ') 



x / e (c 5 (a)-r ) /2)t 2 ^ / e (c,(a)-„' /2)s 2 ^ 



We pick 77,77' > 3c 5 (a). Recall that £ = jy-cN ) / 'dN x / 3 + T 2 . Le t 77 = max(|r|, 3c 5 (a), 1) 
and 77' = max(|r'|,3c 5 (a), 1). It follows from ( |5.32| ) and J5.4C| ) that 

(5.41) \Si\ <C9(/,a), 

z = 1,2 if |r|, |r'| are small. If |r| and |r'| are large, say r, r' ^> 1, then 77 = r and 
77' = t', < r 1 - t < 1, and we get from ( |5 32j ) and ( |5 40fl that 

(5.42) |5i| <c 10 (/,a)e- T3 . 

Inserting these estimates into ( 5.27| ) and using v — u > 1, we obtain an estimate of 
( |5.17 ) of the type we have in the right hand side of (5.2). 
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Consider the expression 

(5.43) E < t ) {x2,XA)fN{x2)fN{x 4 )k l 



X 2 ,Xi£l 



^ab I x 4 x 3 
cd ^ X 2 XX 



In our computations with the kernel K given by ( |5.20 ) we will leave out the complex 
integrations. Thus 



f^ab I x 4 X 3 



K. 



cd 



X 2 Xi 



G* bc (z 2 ,w 2 )^ T G* bd {z 2 ,w 2 )^% 



= G M (2l,M)l)G M (2 2 ,!U2) — — - G ad (z U Wi)G bc {Z2,W 2 )- SI - — . 

z 1 z 2 Zi z 2 

We are led to the symmetrized expression 

(5.44) - E 4>{ x 2,x i )f N {x 2 )f N {x i ) 



x 2 ,x i £l 



Perform the ^-summation first and use the formula (5.25). The parts containing 
Ri and R 2 can be estimated in the same way as above. We will only discuss the 
contraction and finite-difference parts. The contraction part of (5.44) is 



(5.45) 



\ E /»fe) 2 ?[(i+s(^i)r + (i+swn 



x 2 &L 



and hence the contraction part of (5.43) is 

1 \ w X2 w Xl 

7; 22 fN( x 2) 2 {G* a+u _^ c (z 1 ,W 1 )G* bd (z 2 ,W 2 )^—^- 



X2&I 



Ga,c+v-u{ z li w l)G b d( z 2,W 2 )^-^- 



G* a+u - v4 {zi,w 1 )Gl c (z2,w 2 )^4 2 --^- 

z 1 z 2 

G a d ( z i ' w i )Gt c+v - u (z 2 , w 2 ) — } . 



Here we have also used (5.22). Performing the complex integrations we obtain 
(5.46) 

^ E fN( x 2) 2 {K a +u-v,c{ X 2, x 2)K bd (x 3 ,X 1 ) + K a ,c+v~u{ x 2, x 2)K bd (x 3 ,X 1 ) 



X 2 £t 



k a+u - v ^(x 2 ,xi)K bc (xz,x 2 ) - K a . d (x 2l x 1 )K b ^ c+v ^ u (x 3 ,x 2 ). 



The finite difference part of (5.43) is 



(5.47) 



E A/ Ar (x 2 ) 2 {G: c (zi,u; 1 )G^(z 2 , U ; 2 )^%(-+ W i) 



2(1 -a) 1 



21 



%2 %1 1 

G ad (^i, wi)G bc (2 2 , w 2 )- FF - F -(— + w 2 )>. 



-21 
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The double <^-term in (5J5), in ( 5.11 ) and ( 5.14 ) combined give 
(5.48) 



12 53 (j>{x 2 ,x i )cj){x 1 , x 3 )K, 



Xi x 3 
x 2 X\ 



12 ^ ^(11,0:3) 



K v 



3E 



x\ x 2 
Xl x 2 



x 2 x 3 

X 2 Xl 

- 2kr, 



K vu 



Xl X 2 
Xl X 2 



fN(xi)fN(X2)f N (x 3 )f N (Xi) 

fN(xi)fN(x 2 ) 2 fN(x 3 ) 

In(xi) 2 /n(x 2 ) 2 



X 2 X 3 
X 2 Xi 



K: 



Xl x 2 
Xl x 2 



Ai + A 2 +A 3 . 



Consider the 24-summation in A\. The contraction part is, by (5.46), 



fN(xi)fN(x 2 ) 2 fN (x 3 )(j>(xi, X 2 ) 



x 2 x 3 

X 2 Xl 



K v 



x 2 x 3 
x 2 Xl 



which is exactly ~^A 2 . Hence what remains of A 2 is 
(5.49) 



-6^^1,13) 

X 

We have 

E lN{xi)f N {x 3 )K,, 



x 2 x 3 
uu 1 x 2 Xl 



Kt 



x 2 x 3 

X 2 Xl 



fN(xi)f N (x 2 f f N (x 3 ). 



X 2 X 3 
X 2 X 



3 )= E M^)f N (x4)KZ 



xi ,x 4 6 



X4 Xi 
X 2 Xi 



We can now apply ( 5.46| ) to compute the contraction part of the first half of (5.49) 
and get 



(5.50) 



-« E 

X\,X2&t 



K" 



Xl x 2 
Xl x 2 



K uv , 
uv 1 



Xl X 2 
Xl X 2 



fN(xi) 2 f N (x 2 ) 2 . 



Similarly the second half of (5.49) has the contraction part 



(5.51) -6 ^ 



X lt X 2 & 



K!: 



Xl x 2 
Xl x 2 



Xl x 2 
Xi x 2 



]n{xi) 2 ]n{x 2 ) 2 . 



Since the contraction part of ( 5.49| ) equals (5.50) plus ( |5.5l| ) we see that this exactly 
cancels A 3 . It remains to consider the finite difference parts. 
From Ai we get a finite difference part 



(5.52) 
60 



{I -a) 



■E 



G*v u ( Z l . W l ) G 'vu ( Z 2 , W 2 ) -Jjr -Jjp ( — + Wi ) 
Z-^ Z\ 



%2 Xl -1 

-G vu {zi,wi)G vu (z 2 ,w 2 ) — — (—+w 2 ) 

1 ^2 1 



A/at (x 2 ) 2 <K£i , z 3 ) /at (£1 ) /W (a; 3 ) . 
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We also need the finite difference part of ( 5.49 ). These finite difference p arts s hould 
be cancelled by the contraction part of (5.52). The contraction part of J5.49D is 



(5.53) 



3a 



(1-a)' 



E 



w X2 w Xl 1 
z x z 2 Z\ 



- G* vu (zi , wi ) G* uu (z 2 , w 2 ) ( — + w 2 ) 

Zi z 2 Zl 



Af N (x 2 ) 2 4>(xi,x 3 )f N (xiY 



3 a 



(1-a): 



E 



w X2 w Xl 1 

Gt u (zi,W 1 )G* vv (z 2 ,W 2 )^-^r T ( hWl) 

z x z 2 Zl 



X2 X-] i 
f jj ^ f f J L j 

-Gl v (zi,w 1 )G* m (z2,w 2 )- w -^ r { \-w 2 ) 

z x z 2 Zl 

In ( |5.52| ) we have first 

E ~^-^- ( l } ( X T-^ X 3)fN{xi)fN(x3)Af N (x 2 ) 2 

Zi z 2 



Af N (x 2 ) 2 (j)(x 1 ,x 3 )fN(xiy 



^E^!) 2 $(e 



12 „,, x & 



II), W 



w 4>(x 2 ,X 4 )f N (x 2 )f N (x i ) 



This gives the contraction part 

,,Xl „,,X2 



( 5 - 54 ) \ E WfNix^Afxix^iil + g^r-v + il+giw^r-v } 



Z-y Z 2 



From the other half of fl5.52| ) we get similarly the contraction part 

( 5 - 55 ) \ E ^%fN{x 2 ) 2 Af N {x 1 )\l + g{z 2 )r- v + {l+g{w 1 )r^ } 



By (5.54) and (5.55) the contraction part of (5-52) is 

x [G* u (zi, wi)G;„(z 2 ,w2) + G; u (zi,imi)G* 1) (z2,w 2 )] 
3a ^ . , , > 2j - / \2/ 1 

- (1^)2 E ## A ^) 2 ^) 2 ^ + «*) 

x [G* u (zi, wi)G;„(2:2,w2) + G^(2;i, wi)G;„(z 2 , W2)], 

which exactly canc els ( 5.53 ). The finite difference part of ( 5.5 2j ) is handled in the 
same way as (5.28). 

We will end with some brief comments about remaining estimates. By ( 5.20| ) we 
have for example 



(5.56) AK uv {x,y) 



1 



(2wi) 2 
1 



dz / dw 

lr 2 Ay, 



1 — OLZ 



1 — a/z\ N / 1 — aw \ " w y 



1 — a/w J z 



N 



-[(1 + g(w)r~ u - 1](1 + ,gWr(l + g(w)) u . 



w(z — w) 

Here we can expand (1 + g(w)) v ~ u — 1 as in ( |5.24| ) and then estimate in the same 
way as we did for the R\— and R 2 — terms above. In this way we will see that the 
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Aif 2 -terms give contributions of the right type. We get a similar integral expression 
for AK UV + AK VU — AK VV involving 

(1 + g(z)y (1 + 5 M)-"[(1 + g{z)T- v - 1][(1 + g(w)y- u 1] 

and we proceed similarly. 

□ 



5.2. Weak convergence. Consider iJjy(/, tj) as defined by ( |5T| ). The next lemma 
is a standard consequence of lemma 



5.1 



Lemma 5.2. Under the same assumptions as in lemma 5.1, 

(5.57) P[ max \H N (f 1 tj) — H N (f, t u )\ > A] < C(/, a)A~ 4 e~' tl1 ' 3 \t u — t v \ 2 . 

j=u,...,v 

Proof. Let 

r)i = H N (f,t u+i ) - H N (f,t u +i-i), 

T m = YT= l ^ and T o = 0, so that - Tj = ffjv(/, ij) - H N (f, U). It follows from 
( |5.2| ) and Chebyshev's inequality that 

P[|T -Tj\ > A] < Ca^A^e-l'"! 3 I £ U/ 

for u < i < j < v, where it£ = — tg^i +u . This implies (|5.57|) according to 
theorem 12.2 in ||. □ 

Fix I > and consider rescaled top height curve iTv,o(*) for \t\ < L and its 
modulus of continuity, < S < 1 , 

w N (5) = sup \H Nfi (t) - H Nfi (s)\, 

\t\,\s\<T,\s-t\<6 

Lemma 5.3. Let wn be defined as above. Given e, A > there is a 5 > and an 

integer Nq such that 

P[w N {6) > A] < e 

ifN>N . 

Together with the convergence of the finite dimensional distributions, theorem 



1.1, this proves theorem 1.2, We turn now to the proof of lemma 5.3. 

Proof. Assume that (5 _1 ,T € Z and divide the interval [— T, T] into 2m parts of 
length T/m = 5. Write r 3 - = [jS[cN 2 / 3 ]], c = 2a/(l - a), so that t Tj w jS. 

Claim 5.4. Let L = T[c7V 2/3 ] and B M is the subset of our probability space where 
maxui<£ \Hpf (tj)\ < M . Then, given e > 0, we can choose M so that 

(5.58) P[B C M ] < e. 

We will prove this claim below. We will also need 

Claim 5.5. For any A > there is a constant C(M) that depends on M but not 
on A such that 

(5.59) P[ max \H N , (U) - H N . {t r] )\ > \B M ] < 
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We will return to the proofs. The proof of both claims are based on choosing 
appropriate functions / in lemma 5.2 and results about convergence in distribution. 
Assuming the validity of the two claims we can prove lemma 5.3. Set 

(5.60) A 



j - { max |-Hjv,o(s) - H N $(t r .)\ > A/3}, 

t rj <i<t Tj+1 



so that {wm{5) > A} C U\^<Aj. Choose M so large that F[-B^, f ] < e, which is 
possible by claim ^4|. Hence 

(5.61) P[uijv((5) > A] < e + ^ u B m\- 

\j\<m 

Now, if the inequality in (5.60) holds then 

max \H Nfi (ti) - H Nfi (t r ,)\ > A/9. 

Tj <i<r j + 1 



Consequently, using ( 5.59| ) and ( 5.61 



P[w N (S) > A] < e + (2m + 1)C{M )\- 1 8 2 < e + 2TC(M)\- 1 5. 

Choose S so that i5 < e\/2TC(M). Lemma 5.3 is proved 

Consider now claim 5.4. Pick a C°° function q such that < g < 1 and 



(5.62) 



1 if x > 
if a; < — 1 



and let Qm(x) = g{x — M). It is not hard to see that if we take / = gu in ( p.2\ j 
the C(f, a) can be taken to be independent of M (only sup-norms of / and its 
derivatives enter). If ifjv(<7M ,t Tj ) > 1/4, then HN,o{t rj ) > M — 1 and using the 
convergence in distribution to F2 we see that we can choose M so large that 

P[H N (g M) t r} ) > 1/4] <e 2 

for < m and all sufficiently large N. Let uj denote a point in our probability 
space. Now, 

F[max H N (g M ,tj) > 1/2] = F[H N {g M ,t j{u) ) > 1/2] 

\j\<L 

m 

= V[H N {g M ,t j(u) ) > 1/2, < t j(uj) < t rj ] 

j ——m-\-l 
m 

< ]T nH N (g M ,t rjj )>l/4] 

j=—m+X 



E 



j=— m+l 



[ max \H N {g M ,U) - H N (g M ,U ,-_J| > 1/4] 

r-i-\<i<t r . 

J 7 



cs 



by lemma We can now choose 5 = e. 

If Hf{fi(tj) > M, then jffjv(g m , tj-) > 1 and hence mas| 3 |<£ HN${tj) > M, which 
implies max|j|<£ H]si{gMtj) > 1/2. It follows that 

P[max J ff A r (t 7 -) > M] < e. 
b'l<i 



The case maxui<i H^fiitj) < —M, is analogous. This proves claim 5.4 
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To prove claim p.5| we let i(uj) be defined by 

max \H N fl(ti) - H N0 (t r .)\ = |-Hw,o(^(w)) - H Ni0 (t rj )\. 

rj<i<r j + 1 

Let Ij = [jX, (j + 1)A), j = -K, ...,K-l, where M = KX, K e Z+. Take a C°° 
function /, < / < 1, soch that f(x) = if x < -A, /(x) = 1 if < x < A and 
F(a;) = if x > X. Set 



(5.63) 



f j (x) = f(x-Xj). 



Suppose first that -Ha^o^i^)) < HN,o(t rj ) — 2A and that a-> S Bm- Then there is a 
k(uj) such that H^^{t ri ) S /fe( w ), and 

H N , {t i{u)) ) <{k + 1)A - 2A = (fc - 1)A, 

and consequently H N (f k ^, t i(w )) = 0. Since H N (f k ^,t rj ) > 1, we see that 

|-ffjv(/fc(w),£j( w )) - HM{fk(u)jUj)\ > 1- 

Hence, 



(5.64) 



max max \H N (f k ,ti) - H N (f k> t ri )\ > 1. 

fc|<m r 5 -<j<r 3 - + i 



Call this event F. If we instead suppose that -ff/v,o(ii(w)) > HM.o{t rj ) + 2A we can 
proceed similarly and see that (5.64) still holds. Now, 

P[F] < P[U rj < i < rj+1 {|Hj V (/fc,* i ) - H N (f k ,t rj )\ > 1}] 



fe=-/f+i 



A 



by lemma 5.2. 



□ 



5.3. Tran sversal fluc tua tions. In this section we will prove corollary |l.3| , propo- 
sition 1.4 and theorem 1.6. Let T > be fixed and set 

s n( u ) = SU P H N , (t) 

-T<t<u 

S T (u)= sup (A(t)-t 2 ). 

-T<t<u 

We will write S£ for S^(T) and 5 T for S T (T). 
Lemma 5.6. Given e > we can choose T — T(e) so that 
(5.65) F[S% ^S T N ]<e 

for all sufficiently large N . 



Note that together with theorem L2 this proves corollary [1.3] . 

Proof. Let qm be defined as above and set Rj = T + (j — 1)5, j > 1, where S will 
be specified below. It follows from lemma [Tl that for Rj < log N, 



(5.66) 



»[_ sup Iffjv^Af, t) - ffjv(<?M, > 1/2] < Ce^ 2 
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where C is independent of M. Now, 
(5.67) 

P[ sup H N , (t) > M] < P[ sup H N (g M ,t) > 1] 

T<t<R L T<t<R L 

= P[ max sup H N (g M ,t) > 1] 

l <j< L Rj<t<R j + l 
L-l 

< V P[ sup (H N (g M ,t)-H N (g M ,Rj)) + H N (g M ,Rj) > 1] 
~[ Rj<t<R j+1 

L-l L-l 

<VP[ sup (H N (g M ,t)-H N {g M ,R J ))>l/2} + yP[H N (g M ,Rj)>l/2}. 

Claim 5.7. There is a positive constant c\ such that 
(5.68) V[H Nt0 (R) >s}< e -^+R 2 ) 3/ \ 

Proof. Let c = (1 + a)(l - a) -1 ** -1 , d 3 = a(l + a)(l - a)" 3 and ^ = j/cN 2 / 3 as 
before. We have 

1 2a 

F[H N ,o(R) > a] = P[^ T7 ^(/i (2j, 27V - 1) - > *] 



¥[G(N + cN 2/3 R, N - cN 2/3 R) > -^—NsdN 1/3 } 

1 — a 



by (1.4) and the definition of Hj^ q. Recall that the parameter q in the geometric 
distribution = a 2 . By Corollary 2.4 in jl5| we have, for all K > 1 and j > 1, 

(5.69) F[G(tyK],K) > Kt] < e - 2KJ{ - t+1 \ 
where the function J satisfies 

(5.70) J((l + V?7) 2 (l - qy 1 +S)> c[6 3 / 2 

for < 5 < 1; c[ is a positive constant. We take K = N - cN 2 / 3 R, 7 = (AT + 
cN 2 / 3 R)/K and £ = (2a(l - a^iV + sdN^ 3 )/K. Pick <5 so that 

(l + V97) 2 (l-9) _1 +^ = l + i- 



This gives 5 = diV" 2 / 3 (s + i? 2 ) + ©(TV" 1 ) and if we insert this into (|5.70|) , the 



estimate (5.6E) gives us exactly what we want. □ 



If H N (g M ,Rj) > 1/2, then H N . (Rj) > M - 1 and hence 
(5.71) J^nHNigM,^) > 1/2] < I^e-^-i+^+O-D^) 2 )^^ 



< 1 / e-^" 1 ^ 2 ) 372 ^. 



T-l 
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Using ( p.66| ) we find 
L-i 

(5-72) Y, 



sup (H N (g M ,t) - H N (g M , Rj)) > 1/2] 
—[ R 3 <t<R j+1 

L-l 



< J2 Ce- R U 2 < CS 



dx. 



T-l 



Inserting (5.71) and (5.72) into (5.67) gives 



(5.73) P[ sup H Nfi (t) > M] < - 

T<t<R L 

if R L < \ogN. We can take 5=1. 
It follows from (|/n]) that 



3 „ c (Af_l + x 2 ) 3/2 



dx + CS 



T-l 



dx 



T-l 



(5.74) 



sup H N . (t) > M] < P[H Nfi (t u ) > M] 



R L <t 



t u >R L 

t 2 s3/2 



< Y, e-< M+ <y ' < CNe- c ^ N ^ < e/4 

u>cN 2 / a R L 

if N is sufficiently large. We know that P[H Nfl (0) < M] - > F 2 (M) as N oo and 
we ca n choose M so large that the right hand side of (5.73) is < e/4. Together with 
( 5.74 ) this gives (using symmetry), 



(5.75) 



[sup H N . (t) >M]<e. 

\t\>T 



If Hn,o(0) > M and sup| t | >T HN,o(t) > M, then = Sj^- and consequently 
P[S*^ ^ S£] < P[i?jv,o(0) < M] + P[ sup H Ni0 (t) > M] < 2e 



\t\>T 



for all sufficiently large N. 



□ 



We turn now to the transversal fluctuations and the proof of theorem 1.6. Define 

K T N = mi{u > -T ; S&(u) = S%} 
K T = inf{u > —T ; S T (u) = S T }, 

which give the leftmost p oint of maximum in [— T, T] before and after the limit. 



We first prove proposition 1.4 



Proof. (Proposition 1-4)- Note that 



{K N < -T} C { sup H Nfi (t) > H N . {0)}. 

t<-T 



It follows that 



»[K N < —T] < P[H N , {t) > M] + P[-ff/v,o(0) < M} < 2e 



by (5.75) and the discussion proceeding it. Also, {K^ > T} C {S^ ^ Sj^} and we 



can use lemma 5 



□ 
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Proof. (Theorem l.t). It follows from lemma 5.6 that given e > we can choose T 
and Nq so that 

(5.76) P[K N = K%\>l-e 
if N > No. Let h T : C(M) -> K be defined by 

hr(x) = inf{u > — T; sup x(t) = sup ^(i)}, 

-T<t<u —T<t<T 

and let 

Z?h T = { x £ C(M) ; /it is discontinuous at a;} 
It follows from our assumption that P[Z?h T ] = 0, since hr is continuous at x unless 
x has two distinct maximum points. Since Hn,o converges in distribution to X in 
C[-T,T] it follows that 

(5.77) K T N = h T {H N , ) -» h T {X) = K T 
as iV~ — > oo. 

Let X>t be all points of disconituity for x — ► P[AT T < x], T € Z, and I? = 
Ut^iCt- We will prove that 

(5.78) F[K N <x]^ F[K < x] 

as N — > co for all x £ R \ T>, which implies what we want since T> is countable. All 
the results and assumptions that are behind the estimate ( |5.76 ) can also be proved 



for the limiting Airy process and we can assume that N and T e Z + are chosen 
so that also 

(5.79) F[K = K T ] > 1 - e 

if N > N . Let i£l\D. Then, 

F[K N <x}= F[K% <x,Kl = K N ] + F[K N <x.K T N + K N ], 
and similarly for Kjj. Hence, 

\F[K N < x]-F[K% < x]\ < 2e 



if N > Nq. Since x £ M. \ D it follows from ( 5.77 ) that we can ch oose N\ so 
that \F[K% < x] - F[K T < x}\ < e if N > N^ It follows from (|j|) that 
\F[K < x] — P[K T < x]\ < e. Combining the estimates we see that 

\F[K N < x]-F[K < x]\ < 4e 



if N > max(A^ , A^i), which proves (|5.78| ). □ 



Acknowledgement: I thank Peter Forrester for drawing my attention a few years 
ago to the relation between the exponents occuring in |l3| and p0| . 
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